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PAPERS 


DIVERSION TUNNEL AND POWER CONDUIT 
NANTAHALA HYDROELECTRIC 
DEVELOPMENT 


The salient features the diversion tunnel and power conduit the 
Nantahala Hydroelectric Development, North Carolina, constructed 1940, 
1941, and 1942, are described this paper. The tunnel was used carry the 
flow the river during construction, and the finished project provides 
emergency release from the reservoir. self-removing temporary closure 
gate was used. Permanent control attained narrow and long fixed-wheel 
gates operating rigid arrangement guides and liner. energy dis- 
sipator prevents damage the unlined tunnel below the control gates. 

The power conduit carries water from the reservoir the powerhouse. 
There are lined and unlined sections pressure tunnel, pipe line 
diameter crossing Dick’s Creek, fully welded penstock carrying static head 
slightly more than 1,000 ft, surge tank with restricted orifice, rock 
catcher, and intake. The tunnel was located avoid water loss 
seepage. 


GENERAL DESCRIPTION PLANT 


The Nantahala Hydroelectric Plant situated the Great Smoky Moun- 
tains Macon County, North Carolina, about miles west Franklin, the 
county seat, and miles south-southeast Knoxville, Tenn. was designed 
and built 1940 1942, although preliminary designs had been made years 
before. Approximately mile tunnel had been completed 1931. The 
development takes its name from the small river which located. The 
valley narrow and deep, the Cherokees christened meaning 
the the noonday sun.” 

comments are invited for immediate publication; insure publication the last 
discussion should submitted May 1950. 
1Chf. Engr. and Vice-Pres., International Eng. Co., Inc., San Francisco, Calif. 
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Fig. shows the general layout the development. rock and earth-fill 
dam creates the reservoir. order build the dam, was necessary 
drive diversion tunnel. The power conduit extends from the intake the 
reservoir the powerhouse and has total length 29,499.02 ft, 5.58 
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miles, horizontal measurement. Full reservoir water surface El. 2890.0, 
and normal tailwater, El. 1882.0, that the maximum static head the 
plant 1,008.0 ft, The powerhouse contains single Francis-type reaction 
turbine. The dani? and the turbine have been described elsewhere. 


GEOLOGY 


The Great Smoky Mountains form part the Appalachian system and are 
relatively old. the area the Nantahala development, the rock nearly 
all arkose, variety sandstone containing much feldspar. Feldspars are 
silicates aluminum with potassium, sodium, and calcium present varying 
amounts, and the sands were presumably derived from the weathering feld- 


Nantahala Earth-Faced Rock-Fill Dam,” Growdon, Proceedings, ASCE (publication 
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spathic rocks such granite. The rock generally sound, very abrasive, 
and blocky. When work was begun, faults were known the area. 

The overburden red clay, varying widely depth, which stands well 
its natural condition. Previous experience another project similar 
material had demonstrated slides need expected reservoir slopes. 


TUNNEL 


The drainage area above the dam site about miles, and the average 
annual rainfall, about in. The drainage area heavily forested, but 
rugged enough that runoff quite rapid. The average annual flow the 
Nantahala 331 per sec, and the low flow perhaps per sec; but 
flood 60,000 per sec may expected, although none greater than 15,000 
per sec has ever been known occur. was decided provide for the 


Power tunnel 


and access 


Flow—> 


(b) SECTION B-B (c) SECTION A-A 


(d) SECTION ALONG CENTERLINE 


diversion maximum 7,500 per sec. Fig. shows the general 
arrangement the diversion tunnel. 

Four main factors influence the discharge capacity such tunnel: (1) 
The tunnel length; (2) the tunnel area; (3) the head available force water 
through the tunnel; and (4) the roughness the tunnel walls. 
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The tunnel length determined alinement and grade, and these depend 
the terrain. diversion tunnel should low possible reduce the 
height the upstream cofferdam. The entrance must far enough up- 
stream that ample room available for cofferdam construction, and that 
currents will not disturb the cofferdam. The exit must downstream 
sufficient distance prevent the discharge from interfering with the work 
from damaging any structure. possible, the axis the tunnel, both 
entrance and exit, should right angles the surface contours. There 
should also enough rock cover that driving the abutment the dam 
not weakened. 

The remaining three factors are interdependent. The head available 
depends the height which the upstream cofferdam built, and must equal 
the sum velocity head and losses. Relative economics various combina- 
tions cofferdam height and tunnel diameter must considered. The 
Manning formula was used determine losses; thus, 


n)? 


which the friction head loss, feet; the length conduit, feet; 
the quantity water, cubic feet per second; the Kutter roughness 
coefficient; the cross-sectional area conduit, square feet; and the 
hydraulic radius conduit, feet. 

was decided that, after construction, the tunnel would utilized 
emergency outlet, and means providing additional flood discharge 
capacity. Therefore, was necessary design permanent control works. 
Two alternatives were considered—a control tower the reservoir the 
tunnel entrance, and shaft about line with the axis the dam extending 
upward from the the rock surface. The shaft was chosen (see Fig. 
being somewhat more economical and safer. control tower requires forms 
and reinforcing steel, and these may balanced against the cost sinking the 
shaft. Although this area not ordinarily considered subject earth- 
quakes, there has been one disastrous earthquake Charleston, C., and 
major earthquake western Tennessee (resulting the formation Reelfoot 
Lake). Natahala about midway between these places, and control tower 
would more likely damaged earthquake than shaft. 

was obvious that since the control would have capable operating 
under full reservoir head, the control gates and their setting would form 
constriction that would reduce the carrying capacity the tunnel during 
diversion. Therefore, was decided that the gates would installed after 
diversion, while the reservoir was filling, and that, otherwise, much work 
possible should completed before diversion was started. Accordingly— 
the tunnel was driven, the concrete portals were constructed, the tunnel was 
lined upstream from the control gates, and the control shaft was sunk and 
lined before diversion. 

interesting arrangement was made for closure the tunnel after diver- 
sion begin filling the reservoir. closure gate was required which would 
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resist the head while the permanent control gates the base the shaft were 
being installed, and which could easily removed thereafter. The opening 
the intake portal was made high and narrow (10 ft) accommodate 
sliding gate. Fig. shows the design this structure, and Fig. shows the 
closure gate. The latter had skin plate both faces; was fact caisson. 
was determined that not more than 150 head could accumulate while the 
permanent gates were being placed, and the gate was designed accordingly, 
using unit stress 24,000 per was placed portable crane, 
and water pressure held firmly place long the tunnel behind was 
empty. The concrete portal was provided with by-pass valve for filling the 
tunnel after the permanent gates were installed and closed, and with 3-in. 
pipe the roof allow the trapped air escape. The theory was that, when 
the pressures the two sides the gate balanced, would simply float 
the surface. The scheme was successful, but the gate, fabricated, did not 
have strictly planar sealing surface, and permanent gate installation was 
somewhat hampered leakage into the tunnel. 


Closure gate slots 


Filling 


Ei. 2670.0 


SECTION H-H 


Fixed-wheel gates, with cable hoists, were chosen being the simplest and 
most economical control. There are two gates, series; the downstream gate 
considered the operating gate, and the upstream one, the emergency gate, 
although the two are identical. The first installation step was place the 
gate guide and liner assembly, shown Fig. This was very rigid arrange- 
ment, designed four pieces for easy shipment and handling. However, 
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was completely shop assembled and very carefully checked. The field joints 
were bolted, was feared that any welding might disturb the exact aline- 
ment roller paths and seals. The assembly comprised bell mouth, the 
pocket for the upstream gate, tapered section, the pocket for the downstream 
gate, and another tapered section. Air vents were provided downstream from 
each gate, since partial gate openings considerable air required prevent 
irregularities flow. Seal and roller paths are stainless steel. 

Fig. shows the details concrete the shaft. The air vents are extended 
the top the shaft, the downstream vent being made 3.5 diameter 
that may used asa manhole. The vents terminate air passage under 
the gatehouse floor, and open the outer air, order that partial vacuum 
can ever created within the gatehouse. Wood water stops were provided 
all horizontal lift joints. 

The fixed-wheel control gates, shown Fig. were made long and narrow. 
There are two main advantages this shape—the span between wheels 
short, and the gate has less tendency chatter. The gates were filled solidly 
with concrete, and will close their own weight under full head. The rubber 
seals were new design and, when the gates were first closed, these seals were 
absolutely watertight. However, has been found that when the gates are 
operated under full head, the flow over the top the gate and around the sides 
above the ends the seal plates damages the rubber seals. Apparently, the 
seal plates should have been extended upward far enough that flow could 
pass the side seals until the gate was fully opened. The top seal should have 
been protected suitable baffle plates prevent water under high pressure 
and velocity from passing the seal. 

The gatehouse the head the shaft has two stories. The upper story 
carries the gate hoists, and the lower story high enough that completely 
hoisted gate will hang clear the floor and therefore entirely accessible. 
Each gate lifted two wire ropes in. diameter. Hoists are electrically 
driven, and permanent power line the dam was contemplated, prime 
mover was provided the shape gasoline-driven generator set. 
this mountainous country, subject intense thunderstorms, transmission 
line likely give trouble, and the generator set would have been required 
stand-by source power any event. 

was known that only infrequent use would made this tunnel for 
emergency discharge. However, was feared that jet water issuing 
from the gate under 213.5-ft head would tear rock out the walls the 
unlined tunnel. Experiments were made, and energy dissipator was de- 
signed destroy.the energy this jet. Details are shown Fig. Just 
downstream from the gate steel-armored deflector split the jet and throw 
the halves upward and outward. This deflector cylindrical.chamber, 
the lower end which baffle having circular hole diameter. Down- 
stream second cylindrical chamber, baffled turn, but the hole 11.5 
diameter. third chamber follows, the hole this time being diameter. 
The energy dissipator has had one severe test. was operated continuously 
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for more than days. Water issued smoothly and quitely, and later inspec- 
tion showed evidence attack any the concrete walls baffles. 
The rock was pressure grouted from the shaft, and vertical fan-shaped 


pattern extending downward from the area where the permanent control gates 
were later set. 


station 10+82.16 
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dam was located the most economical site, considering 
all factors, such cost, reservoir capacity, etc. Its foundation about 
El. 2670.0 and full reservoir level El. 2890.0. For approximately miles 
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downstream from the dam site, measured along its course, the river flows 
very deep, narrow, and crooked valley, falling this distance El. 1882.0. 
then flows out into wider valley, and its gradient flattens. The mouth 
the narrow valley about 5.5 miles from the dam site, the crow flies. Pre- 
liminary investigation disclosed that the head between the dam site and the 
valley mouth could developed placing the powerhouse the mouth, and 
providing pressure conduit from the dam. The pressure conduit had 
either steel pipe line tunnel, and was easily demonstrated that should 
mostly tunnel. pipe line would have been very long, and difficult and 
expensive construct and maintain. 

Beginning the face the concrete intake, the power conduit comprises: 


Item Description Length, feet 
tunnel 
steel pipe line, across the valley Dick’s Creek, tribu- 
tary the Nantahala (length between portals) 
steel penstock, terminating the upstream face the 
butterfly valve above the turbine 1,247.23 


Total length 29,449.02 


(There surge tank 1,847.23 from the face the butterfly valve.) 


The the Nantahala development was started for the first time 
1930 was decided that the character the rock would permit the use 
unlined pressure tunnels, and general this proved true. was realized 
that the rock was not absolutely watertight, and was obviously more desirable 
have water leak into the tunnels than leak out them. The tunnels were 
therefore located that, any point unlined sections, the static head 
them less than the vertical distance the water table. loss water 
has been observed. The area around the tunnels has been carefully examined 
for any augmentation springs brooks, but nothing the kind has been 
found, and believed that instead losing water the tunnel, actually 
little gained. Alinement and grade have been shown Fig. 

There sharp irregularity the tunnel grade Station 102 92.94. 
Originally, was planned drive from Portal No. adverse grade 
0.25% far Station 119 02.9, and from Portal No. ascending 
grade intersect this station. well hole was sunk the planned high 
point provide air vent. the job progressed, was found 
advisable stop driving from Portal No. Station 102 92.94, and com- 
plete the tunnel from Portal No. final grade was adopted shown. 
Since water velocities are low, the maximum being only 3.66 per sec, the 
irregularity has negligible effect the hydraulic properties the tunnel, 
but does make little inconvenient inspect. 

Four different cross sections tunnel were used, shown Fig.9. Work 
force account was started Portal No. 1930, using section, 
ft, with flat invert and semicircular roof. was found that shot 
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(a) CROSS SECTION (b) CROSS SECTION 
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flat invert and semicircular roof was tried, and was found actually 
cheaper per linear foot than the smaller section. This section was used until 
work was discontinued because the economic depression 1931. When 
work was resumed 1940, the contractor elected drive section, 
ft, with 13-ft flat invert and semicircular roof. The fourth type cross 
section was 12-ft-diameter circle, used where the tunnels were steel lined. 

Tunnel headings were driven full face, with pilot tunnel, except the 
one point where difficult ground was encountered. Station 100 
apparent fault was encountered (extending Station 100 86) which was 
filled with assortment sand, calcite, talc, and disintegrating arkose. 
flow about 150 gal per min water developed. was necessary use 
top heading and bench through this zone, and afterward heavily reinforced 
concrete lining was placed for total length 215 centering the bad spot. 

Although other bad driving was encountered, was necessary place 
plain concrete linings for short lengths four other places where narrow 
zones poor rock were found. 

has been previously stated that the guiding principle design was the 
desirability having water tend seep into the conduit, instead out it. 
This tendency reversed near all portals except the intake. prevent 
water losses, the tunnels were therefore steel lined for varying distances back 
from Portals Nos. and Portals Nos. and the linings were ex- 
tended back points where the ground surfaces directly above are elevations 
slightly higher than full reservoir level. Portal No. the ground surface 
directly over the end the steel lining about 400 above full reservoir level; 
but the mountainside the right the tunnel steep that the minimum 
distance the surface much less. Another reason for extending the lining 
well back this portal that the rock here generally inferior that else- 
where the project. 

The Rock reaction turbine has speed 450 rpm, and any 
solid object passing through with the water would probably cause quite little 
damage the buckets. Also, any sand would damage the turbine seals. 
necessary, therefore, take all practicable precautions the end that 
solids, large small, are carried down the power conduit. 

When all work the tunnel had been completed, the walls and crown were 
inspected carefully and loose material was szaled and, finally, they were hosed 
down. The invert was then very carefully cleaned. The entire conduit was 
finally inspected carefully, from end end, make certain that nothing 
loose had been left it. Then, the last step, the conduit was flushed out 
“cracking” the intake gate, described subsequently. All 
cautions were designed take care any solids debris resulting from 
construction. 

There remained provision made for catching any solids originating 
the tunnel any way entering during operation. previously men- 
tioned, normal water velocity low, and any material moved which not 
buoyant must washed along the invert; but the invert fairly rough, and the 
thousands pockets will ordinarily catch and hold such material. How- 
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ever, whenever the tunnel has been emptied and refilled, there will 
tendency for any loose material move along the invert. Fig. shows the 
rock catcher built intercept such possible loose material just the entrance 
the steel lining within Portal No. The tunnel was widened one side, 
and its invert made lower for length 100 ft. The invert was smoothed 
with concrete, and concrete walls apart were built the level 
the invert the steel lining. Wood baffles apart were placed from wall 
wall right angles the tunnel axis, forming series pockets dead water, 
over which any solid material would have pass enter the steel-lined section 
and the penstock. The baffles were made planks, in. in., with keeper 
bars prevent their floating out the retaining slots. more than 
years operation, apparently solid material has ever reached the turbine, 
and almost material has accumulated the rock catcher. 


End steel lining 


Each baffle 


walkaway side only 


SECTIONAL ELEVATION A-A 
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SECTION C-C 


The upstream entrance the power tunnel the side hill, 
about above the bottom the original river channel. The side hill here 
covered with soil and decayed rock. Previous experience the Santeetlah 
Reservoir North Carolina had demonstrated that this material would not 
slide after submergence, and was therefore decided that should left 
matter economy, and disturbed little possible. The control 
scheme generally similar that the diversion tunnel, and was chosen for 
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the same reasons. The control shaft located about 240 from the entrance. 
became necessary, therefore, design intake structure the entrance. 
disturb the overburden little possible, this concrete structure was 
designed that when completed only long narrow slot would show the slope 
the mountain. The structure two distinct pieces—the one nearest the 
portal being founded solid rock and the outer one the talus slope. De- 
tails are shown Fig. construction adopted permitted bracing the 
walls against each other, and heavy retaining walls were obviated. The tunnel 
lined for distance from the parapet face, point where the 
crown sound rock. 


Center line intake 


(a) PLAN 


rash rack support beams 


Gage well piping 
3in. wrought iron pipe 


Broken stone Construction joint 


(b) SECTION 


INTAKE 


The sloping intake opening covered with steel trash rack, sections, 
laid loosely steel and concrete beams. The racks were designed for 


differential head ft, and the supporting beams for differential head 


ft. The gross area about 532 for contemplated flow 600 
per sec, and the net area, about 70% great. Since the racks can often 
reached during low water seasons, devices for raking hoisting were 
provided. 

The Intake Control.—General arrangement the intake shown Fig. 12. 
Two fixed-wheel gates, series, are placed the foot shaft intersecting 
the tunnel. The gates themselves are very similar the control gates in- 
stalled the diversion tunnel. The guides were not rigidly made, however, 
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and were slightly warped somehow the installation, and was necessary 
very expensive field grinding job the seal and roller paths. The rubber 
seals, similar those the diversion control gates, have proved entirely 
satisfactory here. 


Roller seat guides 


SECTION C-C SECTION B-B 
ATTACHMENT INTAKE GATES 


The concrete wall between the gate wells was made too thin, and therefore 
leaked considerably. has been necessary repair it. 

unusual feature the method the gates—that is, opening 
them only little. When the power conduit has been drained for inspection 
repair, desirable refill slowly. This may done small by-pass 
gate regular gate. The danger using the latter method 
that not foolproof; the gate may opened too wide, and the conduit 
filled too rapidly. The arrangement adopted guards against this possibility. 

Each the gates hoisted two ropes, and will close under its 
own weight. The hoist provided strong enough lift the weight the 
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gate, but not strong enough overcome gate weight plus friction under pressure, 
and the hoist ropes were fixedly connected the gate, could not opened. 
The hoisting force multiplied that sufficient for con- 
necting the hoist ropes pair levers having limited range motion. 
The levers are hinged the gate, but their are the gate guides, 
and the hoisting force multiplied the ratio distances from fulcrum 
hoist rope and hinge, respectively (see Fig. 13). This ratio 5.67. When 


Pinion 


rake 


(b) ELEVATION 


Gate 


hoisting tension put the wire ropes, the long ends the levers move 
through vertical distance in. before they hit rigid stop the gate, and, 
doing, lift the gate 2.47 in. Once the levers hit the stop, the 
hoist ropes effect become fixedly attached the gate, and the hoist not 
strong enough lift farther until the pressures both sides the gate are 
the same. 

Another noteworthy feature simple and foolproof brake the hoist for 
lowering gates. automatic operation, and made standard stock 
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parts. ordinary air blower direct connected the motor shaft, 
such fashion that the blower runs backward when hoisting. Fig. shows 
the arrangement, and its operation explained follows: Assume that the 
gate hoisted, the hoist held stationary solenoid brake, and that power 
available. Release the brake, and the hoist starts overhaul, the gate 
descending. the gate picks speed, the hoist runs faster, and the blower 
begins absorb power. The gate accelerates until develops power the 
same rate the blower absorbs it, and thereafter the speed descent con- 
stant. Since the power that blower requires varies the cube its speed, 
easy design that the speed free descent always less than three 
times the hoisting speed, which usually about all the motor There 
are only two precautions taken—the motor must checked for the 
maximum speed, and the blower discharge and intake must screened that 
they cannot closed accidentally, since closed blower acts the same way 
centrifugal pump, and requires power when does work. 

All parts the hoist are made strong enough stand the stalling torque 
the motor, withstand the maximum pull the ropes. 

illustrate the action the device, and the blower, the 
following data are given: 


Roller and Roller Bearing the formula— 


—the following values applying rollers were: 


Description 

Vertical force, pounds, found necessary over- 

come friction 
Horizontal pressure gate, load wheels, 

pounds 735,000 
Radius wheels, feet 0.645 
Coefficient rolling friction for roller bearings 0.003 
Coefficient axle friction 0.01 
Radius center roller bearings, feet 0.25 


Lift Necessary total weight required the 
gates follows: 


Item Description Pounds 
Weight gate, filled with concrete, 38,700 
Weight wire rope 
Downward hydraulic pressure, top gate 
Roller and bearing friction 
Seal friction, caused mechanical deflection seals, as- 

suming friction coefficient 

Seal friction, caused water pressure seals 
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Dividing 150,660 the lever arm ratio 5.67, the rope pull 26,600 
Thus, the hoist exerts somewhat less force than lifting the dead 
weight the gate. 


Tension Hoist Ropes, Gate Closing Under Pressure, Before 
Device components force considered are: 


Item Description Pounds 


Downward forces exerted gate— 
Downward hydraulic 83,000 


121,700 

Forces resisting downward motion— 
Roller and roller bearing 6,280 
Seal friction, caused mechanical deflection....... 9,130 


Seal friction, caused water pressure 


ee eee eee 


Tension hoist ropes the difference betweenitemsland2 


Tension Hoist Ropes, Gate Closing Under Pressure, Almost Closed, and 
“Cracking” Device component forces for this case are: 


Item 


Description Pounds 
Downward forces exerted gate— 

Pressure overhang head seal 


Forces resisting downward motion (as item C-2) 
Tension hoist ropes the difference between items and 
17,760 lb; divided 5.67 the rope pull 


ee 


The size the hoist ropes was determined item C-3. Two 1}-in. wire 
ropes were used, with total breaking strength Thus, the factor 
safety under emergency closing conditions 2.15. 

The hoist motor must able lift the gate the air, and drive the 
blower backward the same time. Given: Motor speed, 870 rpm; reducer 
230 gear and pinion ratio, 15; second gear and pinion ratio, 
13; and total reduction, thus, 4,140 

The drum speed then 0.21 rpm, and with pitch diameter in., the 
hoisting speed 2.31 per min. The motor must put out, merely hoist the 


gate the air, 3.28 hp, computed follows—P 33,000 0.93 
3.28. This value takes account the fact that the reducer has efficiency 
and each gear and pinion combination, 95%. 


(2) 
000 26,310 
645 
003 3,100 
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650 
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The blower requires 2.0 825 rpm running forward, and one eighth 
much running backward. 870 rpm requires 0.29 running backward, 
3 
computed follows: 0.29. Thus, the motor must put out 
total 3.57 hp. 


The speed which the gate will descend the solenoid brake released 
computed follows: 


Step Description 
Let equal revolutions per minute motor, and therefore blower. 
Speed descent can then expressed, feet per minute, 
or . 
377 
Power generated the falling gate speed descent times tension 
hoist ropes, divided 33,000; thus, 
The power computed step reduced the gear and reducer 
efficiencies 
The power required run the blower 


and 1,344 rpm. This speed conservative, since some 
friction the hoist has been neglected. 


The control house the top the shaft similar that the top the 
diversion control shaft. Power for the hoist furnished small gasoline 
engine-driven generating set. Since the two control houses are only short 
distance apart power line was run between them, that the generating set 
either house may used operate either set gates. Thus, each generating 
set stand-by for the other set. 


preserve reasonable tunnel alinement and grade, was necessary 
come out the open with the power conduit and carry steel pipe line across 
ravine Dick’s Creek, and Fig. shows the arrangement adopted. This 
provided the further advantage additional portals for tunnel driving. Sur- 
face rock this area was poor, and overlaid considerable depth soft 
ground. Sufficient floodway for Dick’s Creek had provided between 
pipe supports. was decided make the pipe diameter the open, 
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and was designed continuous beam over six supports, the longest span 
being ft. Four the six supports were founded concrete caps over 
steel H-beam piling. 

Immediately outside each portal was concentric reducer between 10-ft- 
diameter tunnel linings and the 9-ft-diameter pipe the open. drain valve 
was provided the pipe span over the creek, but because the difference 
diameter the tunnel sections not drain completely, and this objectionable. 
probably would have been better make all this pipe diameter. 

the last stages constructing steel pipe lines and penstocks, access 
their interiors becomes matter first importance. They must thoroughly 
cleaned and painted. After plant has been put operation, such access 
still very important. Repairs are necessary such repainting, cleaning 
out either case, both inconvenient and expensive 
forced pass men and materials through manholes, and easier access 
highly desirable, although manholes should provided for the purpose 
inspection. Dick’s Creek 14-ft length pipe was made removable. The 
reducer the stream portal forms part the removable section. each 
end the section flanged (see Fig. 16) and the flanges drilled for 
stud bolts about in. apart the pitch circle. corresponding flange 
placed the end the 10-ft-diameter pipe upstream. flanges are grooved 
for round rubber packing. the detail the downstream end this 
removable section which makes easily workable. The section flanged 
pipe adjacent only in. long, and terminates expansion joint. When 
the expansion joint loosened, and the flange bolts are removed, the short 
section pipe can slid away from the removable section, that the latter 
has ample room come out. only necessary put blocking under it, 
take few turns cable around it, and rollit was easy make 
this provision; unquestionably, paid for itself during construction, and 
will probably pay for itself several times over the course one repair job. 


The penstock considered extend from point the tunnel where steel 
lining begins the upstream flange the butterfly valve the powerhouse, 
which lies the foot spur jutting out from Rattlesnake Mountain. That 
part the penstock the open runs down along the narrow spine this spur. 
Pertinent data steel pipe given Table distances given being along the 
penstock and not horizontal, and Fig. shows the general arrangement. 

Table shows the location three removable sections. The upper one 
generally similar the one Dick’s Creek previously described, and pro- 
vided for the primary purpose access. The middle one much the same, 
except that not reducer, and was provided essentially for different 
reason. Although may afford access, was essential the scheme adopted 
(described subsequently) for flushing out the tunnel invert. The third and 
lowest removable section adjacent the butterfly valve the powerhouse, 
and quite different from the other two (see Fig. 16). 

was not practicable put expansion joint here, but was necessary 
provide some room for possible removal the butterfly valve. The last 


in. holes at approximately 9 centers for 
lin.¢% x 10 in. bolts with cquare headc 
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section rigidly held penstock was provided with flange, from the valve 
penstock. The removable section flanged match this inclined flange 
one end and the valve flange the other, and thus slightly wedge shaped, 
shown Fig. 16, the wide part the wedge being top. under the 
powerhouse crane, and can lifted straight out. carefully machined 
close fit; will shrunk using dry ice when removal necessary, and 
was originally placed using dry ice. The packing the joints placed from 
the interior the pipe after the removable section place. The flanges 
are recessed shown; the packing placed the recesses and held with 


retaining ring sprung into place and tack welded. This type joint has been 
very satisfactory. 


Wall Center-line Maximum Wall Center-line Maximum 
(inches) (feet) (feet) (inches) (feet) (feet) 


(1) (2) (3) 


(4) 


(1) 


(2) (3) (4) 


0.55 210 13/8 2,353.52 536.48 
1/2 128 2,427.50 462.50 17/16 2,324.00 566.00 
9/16 11/2 2,291.70 598.30 
5/8 19/16 2,262.64 627.36 

19/ 2,236.80 


Fic. 17.—Prorime or Pensrock AND Sur@E TANK 


(b) 8.5 57.5 1,995.90 894.10 


upper end run. Anchor Portal No. vertical and horizontal point intersection (P.I.). 
section. Reducer. Anchor vertical P.I. Vertical P.I. inside powerhouse. 


The middle removable section the penstock was left out when the pen- 
stock was built. After all work the surge tank and power conduit above 
this point (including painting) was completed, and—so far was known—the 
interior had been cleaned scrupulously, the tunnel invert was flushed out after 
“cracking” the intake gate. The tunnel then filled far the apex, and 
since the water velocity upstream from the apex was low, the flushing did 
good this part the tunnel far moving anything heavier than water 


3/4 384 
13/16 
15/16 
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was concerned. was expected, however, that this procedure would float out 
any loose pieces wood overlooked cleaning; but, from the apex down- 
stream, the water flowed turbulent stream down the invert. 

Where the removable section had been left out, temporary funnel section 
was attached the penstock, and the discharge was piped over the side the 
mountain where could damage. Water was allowed run this 
way for days, and the effluent was carefully observed. first, this effluent 
was quite murky because the rock dust carried. Any solid object could 
heard striking the funnel. From these sounds, was obvious that the solid 
particles were all small, with one exception. spite the rigorous inspection, 
one piece wood, in. in., and several feet long, came down. When the 
water became clear, flushing was discontinued, the funnel removed, the re- 
movable section inserted, and the power conduit was ready for service. 

The penstock was anchored carefully four points—at the portal, each 
the two vertical curves, and the lowerend. The anchors were designed 
resist hydraulic thrust, and carry the dead weight pipe. asteep slope, 
the latter becomes appreciable item, since all longitudinal forces are carried 
the anchors. For example, the anchor Station (see Fig. 15) 
resists any longitudinal force because the weight pipe between the expansion 
joints near the portal and about 370 below the anchor. the lower end, 
just above the powerhouse, the pipe was incased concrete. The incasement 
provided continuous anchor about 238 long. Its purpose was insure 
that movement the penstock could possibly take place which might con- 
ceivably disturb the alinement the generating unit. 

Ordinarily the penstock will drained through the turbine, but was 
deemed necessary provide for draining case the turbine cannot run. 
The penstock just inside the powerhouse, and 12-in. pipe, controlled 
gate valve, taken off and connected manifold mounted corner 
the powerhouse. The manifold feeds four nozzles diameter, the 
nozzles are inclined 45° with the horizontal, and point the Nantahala 
Valley the right of, and upstream from, the powerhouse. Each nozzle 
controlled gate valve. When operated, the jets from these nozzles tear 
themselves pieces the air, and the discharged water falls heavy mist 
over area perhaps 200 wide and 600 long. determined that more 
than 6,000 dissipated this way and, since all done air resistance, 
sets considerable windstorm the valley. 


TANK 


The surge tank was placed close possible the powerhouse, consistent 
with not having too tall tank projecting above ground surface. Because 
the surface contours, was not advisable locate directly above the tunnel, 
and was actually slightly more than the left the tunnel axis. The 
riser was inclined with elbow the bottom, shown Fig. 
arrangement permitted the use rock bin the bottom the riser, that 
footings for the riser and surge tank could excavated without interfering 
with work the tunnel. 


tio: 
Iten 
dis 
cal 
a 
tan 
thi 
exc 
anc 


December, 1949 NANTAHALA DEVELOPMENT 1435 


The surge tank 600 from Portal No. and the most important eleva- 
tions are follows: 


Item Location Elevation (feet) 
Center line tunnel this point...................... 2,426.50 
Top riser, and bottom surge tank................. 2,715.00 
Rock surface, after excavating flat 2,840.00 
Full reservoir water surface 


The riser steel lined, with inside diameter ft. The surge tank has 
inside diameter 19.75 ft, and has total height 230 ft, the lower 125 
its height being shaft excavated the mountain. Because the 
distance between the surge tank and the powerhouse, speed and load regulation 
delicate. 


in. plates 


(b) SECTION A-A 


18.—Port Tank RISER 


Considering all conditions, either differential surge tank one with 
restricted orifice was indicated. The differential tank would have required 
bottom El. 2630.0 and, course, internal riser. The other type 
tank required the development restricted orifice with coefficient dis- 
charge for water passing upward about 50% greater than for water passing 
downward. Experiments were made various types orifices and one 
finally was designed with discharge coefficients 0.85 and 0.60. Fig. shows 
this orifice port. Computations indicated that, with this orifice, tank 
could built with its bottom much higher elevation. The savings 


excavation and steel made the orifice type tank much cheaper, and therefore 
was used. 


AND WELDED STEEL AND SurGE TANK 


The steel for all tunnel linings, the Dick’s Creek pipe line, the surge tank 
and its riser, and the penstock, conformed the Specification for Flange Steel, 
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Serial No. A-70-39, American Society for Testing Materials. This steel was 
chosen because its suitability for all welded construction. welding con- 
formed the Code for Unfired Pressure Vessels (1938 edition) the American 
Petroleum Institute and the American Society Mechanical Engineers. 

This code permits working stress 25% the minimum tensile stress 
steel, multiplied four factors, follows: 


Item Material Factor 
For flange steel 
For double-welded butt joints 
For stress relieving 
For radiographing 


The permissible working stress thus the product these four factors and one 
fourth 55,000, and proves 12,667 Actually, 12,600 per 
in. was used, designing for instantaneous stress resulting from the most 
severe combination static pressure, water hammer, and surge. pipe 
wall thicknesses thus determined, in. was added corrosion allowance. 

Two purposes were kept mind designing tunnel linings. These 
were the resistance stress and the prevention leakage. Obviously, the 
point pipe enters tunnel derives added strength from any surrounding 
rock. Also, some point within the tunnel, good rock becomes capable 
taking all the stress. Theoretically, then, the pipe the entrance must 
capable, itself, resisting all bursting pressure, and may progressively 
made thinner until zero thickness required. Practical considerations pre- 
clude theoretical design. one knows exactly how compute the resistance 
that rock offers bursting pressure, and pipe must have certain minimum 
thickness make stiff enough handle. 

Conservative design assumptions were therefore adopted. was assumed 


that all bursting pressure would safely taken rock any point the 


tunnel where upward pressure was balanced the sheer rock in- 
cluded between two diverging lines. These lines originated the ends the 
horizontal pipe diameter, extending upward and outward angle 10° 
from the vertical. Between such points and the portals, linings were designed 
capable themselves resisting all pressures, and thicknesses deter- 
mined the corrosion allowance in. was added. Beyond these points, 
pipe wall thickness was decreased in. per course minimum in. 
Lengths such courses varied from ft. The minimum thickness 
lining was originally intended carried points which the vertical 
distance the. water table would exceed the static head, with the idea that, 
parts the tunnel thus left unlined, water would tend seep rather than 
out, previously stated. However, all three places where steel linings were 
used, they were extended farther than originally intended. Portal No. 
the lining was extended 414 ft; the decision extend was made during con- 
struction and plate was not readily available. Plates 0.41-in. thickness 
and 0.55-in. thickness were used. 

For pipe the open air, the design spans and supports was guided 
the results tests reported Templin and Sturm, Members 
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ASCE, These tests indicated that the greatest stresses pipe 
span occurred supports when was just full water under additional 
head, and that, pipe stiffened adequately supports for this condition, 


in. flange plate 


plate 


wide flange 


concrete foundation 


END VIEW 
TYPICAL SADDLE AND SUPPORT 


may otherwise considered ordinary cylindrical beam. The tests 
also indicated special arrangement supports that would cause minimum 


Stresses. Fig. shows typical details. The pipe has stiffeners only sup- 
Transactions, ASCE, Vol. 98, 1933, 154. 
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ports. Pipe was assumed act continuous beam with break ex- 
pansion joints. Longitudinal stresses were computed be, general, quite 
low. The maximum the penstock was only 4,000 per in. and the 
Dick’s Creek crossing only 6,500 Strain gage measurements made 
the field, with and without water the pipe, indicated stresses satisfactory 
agreement with computed stresses. 

Every precaution was taken insure good welding the shop and the 
field. Pipe was shop assembled sections not exceeding length 
tons weight. Longitudinal welds, which are the most highly stressed, were 
thus generally made the shop (the exception being those the surge tank) 
and were stress relieved and radiographed. Circumferential welds were stress 
relieved. shop welds were machine made double V-butt 
joints. Circumferential shop welds were hand welded double V-butt joints. 

Pipe sections were shipped the job with interior spreaders preserve 
their shape. especially important that sections kept round and true 
the ends, where they butt against adjacent sections for welding. 

Circumferential field joints pipe the open were double V-butt welds. 
Welds all plate in. thick and greater were radiographed. Where the pipe 
shell was in. thick more, was electrically heated 300° before welding. 
radiograph joint, radium was placed the center line the pipe, and 
photographic film was placed outside, and contact with, the joint which 
picture was desired. Lead shields were hung inside the pipe either side 
the radium protect workmen from the dangerous radiations. 

circumferential field joints tunnel linings, backing-up strips single 
V-butt joints were used, they could only welded from the inside the pipe. 

Interior surfaces all pipe steel were given one coat bitumastic solution, 
applied cold. Exterior surfaces were not shop painted. After steel had been 
placed and welded, interior surfaces were given additional cold coat 
bitumastic solution, and final coat bitumastic enamel was applied hot. 
Exterior surfaces against which concrete was placed were left unpainted 
and exterior surfaces exposed the air were given one coat alkyd base red 
lead primer, followed two coats aluminum alkyd vehicle. After the 
development had been operation years, the paint was still “standing 
up” well. 

Steel tunnel lining was placed sections various lengths, the maximum 
being ft. The pipe sections were shop equipped with grout connections 
the top and Those the top were spaced about centers, 
and 2-in.-diameter pipes were inserted make contact with rock. Those 
the bottom were spaced about centers, and were for pipe. After 
concrete had been placed around the steel lining, and had set thoroughly, the 
2-in. pipes the top were grouted, under pressure rising gradually 200 
per in. Quite often the concrete placed under the invert does not make 
good contact with the steel, probably because air bubbles, and the 
pipes were grouted low pressure insure contact. High pressure bulges 
the steel very quickly. 

each the three places where steel tunnel linings ended, concrete 
lining was placed for beyond the steel. radial holes 45° intervals 
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were drilled into rock through this lining, and the holes grouted 500 
When this grout had set, second set eight holes, inter- 
mediately spaced, was drilled and grouted the same pressure. Thus, 
barrier was created against seepage along the tunnel. 

The steel tee connecting the surge tank riser and the tunnel lining was shop 
fabricated, with all welds stress relieved and radiographed. known that, 
under pressure, tee will develop stress concentrations, and was considered 
inadvisable accentuate these field welding. 

The surge tank riser and the surge tank lining were grouted after concrete 
backfill had set insure that hollows existed which would permit the 
lining under variations water level. 

All pipe anchors were designed for safety factor against overturning, 
and were dimensioned that the force tending slide the anchor was not 
more than 0.8 its weight, considering the anchor resisting (by gravity 
alone) the most severe combination forces that could imposed. 
addition, all but one were doweled into rock. One anchor was founded 
soft brownish shale and was not doweled; and the maximum foundation 
pressure was limited, this case, 2.5 tons per ft. 

Surge tank and penstock designs were based the assumption that tur- 
bine discharge 700 per sec would shut off sec. The effect 
water hammer was computed the arithmetic integration method introduced 
Norman ASCE. Computations for surge were also made 
method arithmetic integration. Efficiency tests were run the turbine 
August, 1942, using the Gibson method, the action the surge tank being 
observed the same time. Results indicated that the action accordance 
with design assumptions, and may predicted for any operating condition 
with reasonable degree accuracy. 

Observations friction head loss all rates flow indicate that Kutter’s 
has values 0.011 for lined tunnel steel pipe, and 0.0295 for unlined 
tunnel, respectively, determined the Manning formula. 
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NUMERICAL COMPUTATION BUCKLING 
LOADS FINITE DIFFERENCES 


SYNOPSIS 

After introducing the concept central differences function ap- 
proximate its derivatives simple numerical expressions, the writer shows 
how problem involving linear differential equation the homogeneous 
type, with homogeneous boundary conditions, may reduced problem 
involving the solution system simultaneous linear algebraic equations. 
this manner solutions buckling problems may obtained purely 
numerical computations, using procedure successive approximations, which 
greatly enhanced simple method extrapolation. The errors involved 
the procedure are discussed and complete tables are given used 
connection with extrapolation. The procedure then applied various 
cases buckling beams, plates, and shells demonstrate the technique 
its application and solve few problems not appearing the literature. 


INTRODUCTION 


The evaluation buckling loads beams and plates often presents 
analytical problem, whose rigorous solution may become practically 
impossible when the shape the plate irregular when the variation the 
moment inertia the beam plate does not follow simple law. 

For these reasons, approximate methods computation have been devised, 
the best known among them perhaps being Rayleigh’s energy method,? which 
gives upper bounds for the buckling loads. With the help other methods 
giving lower possible squeeze the approximate load between 
given limits. 

comments are invited for immediate publication; insure publication, the last dis- 
cussion should submitted May 1950. 


Prof. Civ. Eng., Columbia Univ., New York, 


Elastic Timoshenko, McGraw-Hill Book Co., Inc., New York, 
Ist Ed., 1936. 


Gegenstiick zum Ritzchen Verfahren,” Trefftz, Proceedings, International Cong. for 
Applied 1926. 


Minimal Problem Theory Weinstein, Journal, London Mathematical 
Vol. 10, 1935. 


des équations des plaques Weinstein, Mémoires des Sciences Mathématiques, 
ol. 88, 1937. 
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Purely numerical methods for the solution structural problems have 
become very popular among engineers, particularly since 1919, because their 
simplicity and because the wider availability calculating machines. 
Hardy Hon. ASCE, and have pioneered this 
field. 

The numerical and graphical computation buckling loads successive 
approximations has been introduced Vianello and Stodola, connec- 
tion with struts, and their method has been greatly improved 
Newmark,® ASCE. Mr. has applied his relaxation procedure 
the same problem. Since 1910, however, has suggested 
finite difference procedure for the determination characteristic values 
differential equations, which was investigated independently 
This procedure seems present very wide range applications well 
great simplicity. the aim this paper give comprehensive presenta- 
tion the finite difference method, together with some extrapolation pro- 
cedures, greatly enhancing its usefulness. 


CENTRAL DIFFERENCES 


Let the values Taylor expansible function f(x) known evenly 
spaced points the z-axis abscissas: 
3h, (Fig. and also, for the time being, the middle point 


the intervals width defined. The “first central difference” 
indicated the operational symbol (first introduced Shep- 
defined 


The “second difference”’ defined the difference the first 
difference and indicated the symbol f(z): 


Continuous Frames Distributing Fixed-End Hardy Cross, Trans- 

actions, ASCE, Vol. 96, 1932. 
Methods Engineering Science,” Southwell, Clarendon Press, Oxford, Eng- 

and, 1940. 
and, 

Procedure for Computing Deflections, Moments, and Buckling Loads,” 
Newmark, Transactions, ASCE, Vol. 108, 1943. 

Approximate Arithmetical Solution Finite Differences Physical Problems Involving 
Differential Equations with Application the Stresses Masonry Richardson, 
Philosophical Transactions, Royal Soc. London, Vol. 210, 1911. 

Bericht: Berechnung von Eigenwerten,” Collatz, Zeitschrift 
fir angewandte Mathematik und Mechanik, Vol. 19, 1939. 
Differences Formulae,” Sheppard, Proceedings, Mathematical Soc., 

ol. 31, 1899, 
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“Higher order differences” are defined the recurrence equation: 


which leads 


and 


The numerical coefficients are the coefficients the binomial 
expansion and even order differences involve the values f(x 
whereas odd order differences require the knowledge 
+th/2). these last values f(x), Mr. introduced 
new operator called the “averager” and defined the equation: 


whose 


and whose nth power— 


—has numerical coefficients identical with the coefficients the binomial 
expansion divided 2". Comparison Eqs. and shows that 


f(x) or, writing this equation purely operational form 


dropping the symbol both its sides, that 


The successive derivatives the function are represented subsequently 
the differential operator 


whose nth 


may easily shown that the three operators and satisfy, formally, 
the fundamental laws algebra and that hence they can used they 
were algebraic quantities, provided that the result the operations inter- 
pretable terms their definition. 
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Some important relationships between the operators and may 
easily obtained from the Taylor expansion f(z): 


which, the function and remembering the power series expan- 
sion e*, may written symbolically 


(9a) 
Changing into Eq. gives 
(9b) 


operational form, 
(10a) 


Expanding sinh into power series, the operator can finally 
represented the following series involving all the derivatives f(x): 


symbolical form, 
Multiplying Eq. 106 Eq. 11, 
D5 hi D7 
and, squaring Eq. 11, 
hs Dé hs D8 


or, 


mul 


and 
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2m 
Thus, the operators and being integer) may expressed, 


terms the derivatives power series involving only the even 
powers 

The inverse relationships, giving the successive derivatives terms 
its differences, are obtained follows: Solving Eq. 10a for and remem- 
bering the power series expansion 


or, making use Eq. eliminate all the even powers 


multiplying Eq. Eq. 14a, 


5 7 
and, squaring Eq. 14a, 
8 


DIFFERENCE OPERATORS CARTESIAN COORDINATES 


Ordinary the approximate solution differential equa- 
tions finite differences often assumed that all the terms the right- 
hand members Eqs. beyond the first are negligible and that, hence, 
the derivatives may approximated 


relations imply that the higher derivatives are not large comparison 
with its lower derivatives and that the higher powers are negligible 


comparison with its lower powers. (The error inherent this approximation 
evaluated Section 5.) 


Similarly, squaring Eq. and again using Eq. 
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The right-hand members Eqs. are called the order central 
differences” and may represented graphically the convenient 
schemes Fig. 

taking into account two terms the series Eqs. and the deriva- 


2.—First 


follows: 


The right-hand members Eqs. are represented graphically Fig. 
Second order difference operators are more accurate than first order 
operators, but involve the values f(x) larger number points. The 
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Nth order difference operators are similarly obtained. Since difference ratios 
approach the corresponding derivatives approaches zero, the accuracy 
the numerical approximations may increased either making smaller 
taking higher order differences. Which procedure more advan- 
tageous depends the type problem hand, but seldom that the 
additional labor required higher order differences warranted practice 
the results obtainable. 

(Expressions for the derivatives terms differences may also obtained 
passing nth order parabolas through the points abscissa and 
ordinate and substituting the derivatives the inter- 
polated parabola for the derivatives f(z). When the interpolation points 
are symmetrically located with respect the formulas thus obtained are 
identical with the central difference operators given Eqs. and 16. When 
the interpolation points are not symmetrically located with respect the 
formulas give the derivatives terms This latter 
type formula conveniently used the boundaries beams and plates, 
but will not considered this paper. Mr. Southwell most the 
central and lateral difference formulas needed the applications, under the 
name 


DIFFERENCES 


Partial first order and second order difference expres- 
sions, Eqs. and 16, may used evaluate the partial derivatives 


constant computing the differences. Thus, using first order differences 
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and calling the spacing points the y-direction, 


and 


whereas the second mixed derivative with respect and 


f(x,y), obtained applying consecutively the operators the right- 
hand members 17: 


4 
and the fourth mixed derivative squaring the operator D,,: 


Similarly, the first order difference operator corresponding the Laplacian 


Using second order differences (Eqs. 16) the expressions for second order 
partial difference operators may similarly obtained. 

The schemes Fig. and Fig. give some the first order and second 
order operators most commonly encountered, showing how partial difference 
operators involve the values the function the corners rectangular 
square grid. Partial operators Cartesian coordinates are well adapted 
the solution problems involving rectangular square domains. 
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(a) SECOND-ORDER 
OPERATOR 


(c) SECOND-ORDER 
OPERATOR 
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DIFFERENCE OPERATORS SKEW COORDINATES 


The solution problems involving rhombic domains (for instance, skew 
plate problems) obtained more conveniently oblique Cartesian, 


Fig. indicates the relationship between the Cartesian coordinates 
point and its skew coordinates u,v: 


and 


(in which the angle skew) from which easy derive the inverse 
relationships: 


and 
tana 


Given function its derivatives with respect and may 


and Using Kq. 


calcul des plaques obliques par méthode des équations aux differences,” Favre, 
Mémoires, Association Internationale Ponts Charpentes, Vol. 1943-1944. 

Bereuter, Publication Laboratoire l’Ecole Polytechnique Fédérale, 1946 
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4 


The expressions for the higher derivatives are obtained repeated ap- 
plication 21: 


The corresponding difference operators are derived substitution 
central difference ratios taken the and the v-direction for the 
derivatives the same directions, according Eqs. 16, depending upon 


(a) FIRST-ORDER SECOND-ORDER 
OPERATOR OPERATOR 


whether first order second order operators are required. Thus, letting 
and the spacing the points the and the v-direction, 
then, with first order differences, 


rh) 


DIFFERENCE OPERATORS TRIANGULAR COORDINATES 


When the problem solved involves domain coverable equi- 
lateral triangular grid (equilateral triangular plate, regular hexagonal plate, 
etc.), convenient use the coordinates” shown Fig. 
which point located coordinates measured along three axes 
and making 60° angles between them. 


022 1 ‘ 02z 2 
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Squaring Eqs. 24, the second derivatives with respect and become 


2 
Adding both 


sides the latter equation, the expression for equitriangular coordi- 
nates, 


derive the corresponding difference operator, substitute Eqs. for 
the derivatives the z-direction, the u-direction, and the v-direction the 
corresponding difference ratios the same directions. the equal 


spacing the corners the grid the three directions, using first order 
differences, 


(b) 


and 

60° 

and 

2 2 2 2. 

can 
corre 
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The graphs Fig. give the expressions first order and second 
order differences and first order differences. 

ASCE, has given the expressions for and V‘z 
triangular which finite difference operators 


FIRST-ORDER 
OPERATOR 


(c) SECOND 
ORDER 
OPERATCR 


(b) FIRST-ORDER 
OPERATOR 


can applied domains coverable grid equal irregular triangles. 
Mr. Jensen’s finite difference operators are first order operators. 


ERROR THE DERIVATIVES DUE THE SUBSTITUTION 
DIFFERENCE RATIOS 


noted Section the solution differential equations finite 
differences the derivatives the unknown function are approximated the 
corresponding difference ratios. 

When first order differences are used, the error e’; the value the first 
derivative immediately found from Eq. 10b: 


Skew Slabs,” Jensen, Bulletin Series No. 332, Univ. Eng. Experiment 
Station, Urbana, 1941. 


1 
r 
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Similarly, the error the second derivative becomes Eq. 11: 


and, more generally, the error the nth derivative the type: 


which the functions are independent 

Eq. shows that the error due first order differences depends the 
even powers only and that, when the Taylor series f(x) converging 
rapidly (small values and rapidly decreasing values the successive 
derivatives f(x)), the series its right-hand member will also converge 
rapidly—the more becomes greater. 

evaluate the error due second order differences, note that Eq. 16a 
the first derivative approximated the difference ratio: 


Substituting this latter equation the expression obtained from Eq. 
and the expression obtained from Eq. 12a, the error becomes 


Similarly, Eqs. 16b, 11, and 126, the error the second derivative 


and, more generally, the error the nth derivative the type: 


‘When the Taylor series rapidly convergent, the series the right- 
hand member Eq. 29a also rapidly convergent, and more than the 
series Eq. 28. 

The derivations this section show that the error due the approximation 
the nth derivative mth order difference ratios the type: 


which the functions are independent 

the basis this result (due Mr. Mr. was 
able prove that the error the solution problem defined linear 
differential equation and linear boundary conditions, due the substitution 
difference ratios for derivatives, also the expressed 
Eq. 


ential 
the 
the 
opera 


which 
depen 
and 
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ERROR THE CHARACTERISTIC VALUES DUE THE 
SUBSTITUTION DIFFERENCE RATIOS FOR DERIVATIVES 


The solution buckling and vibration problems requires the determination 
the values certain parameter for which the linear homogeneous 
differential equation the problem has solutions, satisfying linear homo- 
geneous boundary conditions, and not identically equal zero. 

Calling one such solution and indicating and two linear 


combinations derivatives buckling equation, for instance, may 
simply written 


Next indicate and the difference operators corresponding and 
solving the difference rather than the differential problem; and the 
error that is, the difference: 


With these symbols the approximating difference equation becomes 


Indicate and the errors due substitution finite differences for 


and 


from which 


(34) 


his classical paper finite differences applied the solution differ- 
ential equations, Mr. makes the assumption that the error 
also the proves this assumption generally untrue. 
the whereas general, Therefore, performing 
the division indicated Eq. and assuming, for instance, first order difference 
operators, takes the form: 


which not the type, both because the first term and because the 
dependence the other (This dependence due the fact that 
and are obviously functions the spacing h.) 


4 
; 
' 
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THE EVALUATION BUCKLING LOADS 
FINITE DIFFERENCES 


When the differential equation— 


—of problem known, together with the corresponding boundary 
conditions the characteristic values may determined the 
following method successive approximations, based the use finite 
differences: 

The domain under consideration covered one the grids considered 
Sections and some whose corners will fall the boundary the 
domain (boundary points) and some which will fall inside the boundary 
(body points). The values the boundary points and some neighboring 
points will determined the boundary conditions the problem, directly 
terms the body points, will noted detail from the applications 
Section 10. The values the body points the mesh ---, 
are unknown instead. 

The difference 


—corresponding the differential Eq. satisfied each body point. 
Hence the application Eq. 37a each the body points gives rise 
set linear, algebraic, homogeneous equations the unknowns f;: 


and this system has solutions not identically zero, and only its deter- 


the degree the determinantal Eq. (k-equation), its roots 
give approximate values the first characteristic values, the original 
differential equation, Eq. 36. problems elastic stability the designer will 
mainly interested the lowest buckling load, that is, the lowest char- 
acteristic value, and will hence solve Eq. for its smallest root Once this 
done for value the mesh size obtaining value for the com- 
putations may repeated for new mesh size smaller than obtain 
better approximation reduction the mesh size and correspond- 
ing increase the the body points considered, may evalu- 
ated, theoretically, any degree accuracy. 

Instead, increase the number body points actually increases very 
rapidly the labor involved writing out and solving the determinantal equa- 
tion and becomes practically impossible better the accuracy the solution 
beyond certain limits. For this reason important use practical 
method evaluation for the determinantal equation and simple procedure 
for its solution. 
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Once the determinant has been obtained the schemes Sections 
and the simplest procedure for writing the determinantal equation the 
method condensation.” can easily that the value 


eee eee 


which each element second order determinant with the same leading 
element called the Any other element a;; may used 
pivot, but this case the determinant Eq. must given the sign 

successive applications Eq. 39b Eq. the order the deter- 
minant A(k) reduced and the k-equation thus obtained. this 


evaluation the multiplying factor may dropped each step, since 


the determinant must equated zero. All the elements any row 
column may also divided common factor for the same reason, thus 
simplifying the numerical coefficients the final k-equation. 

When determinant high order, each term the successive 
determinants obtained pivotal condensation high degree polynomial 
and writing out the k-equation may become cumbersome. easier, 
such cases, substitute approximate numerical value A(k) and 
evaluate pivotal condensation the corresponding numerical determinant. 
Two three trials will usually give acceptable result after interpolation. 

The most practical method for the determination the lowest root the 
determinantal equation Newton’s method the method the tangents,” 


EXTRAPOLATION PROCEDURES FOR THE EVALUATION 
BUCKLING LOADS 


shown Section the approximate values approach the true value 
the mesh size approaches zero. The mode this approach depends 


and Matrices,” Aitken, Oliver and Boyd, London, 1939, 45. 


Mathematics for Engineers and Physicists,’’ Sokolnikoff and Sokolnikoff, 
McGraw-Hill Book Co., Inc., New York, Ed., 1934, 32. 


Gin 
—is equal the value the following determinant order 
eee eee eee 
al 
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the particular problem hand—that is, the differential equation and 
the boundary conditions, and not known, offhand. 

the basis extensive experimental investigations the numerical 
evaluation buckling loads finite differences, the following qualitative 
statements may made the behavior 


many buckling problems, for which exact solutions are known, the 
approximations (k) approach the correct value monotonically (Fig. 10(a)) 
—that is, either continuously decreasing continuously increasing. 

some cases the approach not monotonic for large values 
but becomes rapidly monotonic approaches zero (Fig. 


with decreasing amplitude (Fig. 10(c))—that is, approaches alternatively 
from above and from below. 

some cases the approach one sided (from above below), 
but not monotonic (Fig. 10(d)). 

other kinds approach have been encountered thus far. 

choosing proper sequence values for possible all cases 

The error corresponding monotonic sequence found, empirically, 
the is, for monotonic sequence the functions 
Eq. are practically independent and practically zero. The 


ce 
are 
sub 


ases 


ally, 
n(h) 
The 
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error e’, due first order differences, may thus written 


This last result makes possible evaluate the characteristic value 
extrapolation procedures, first suggested Mr. which enhance 
the powerfulness the method considerable extent. 
Eq. 40a, when sufficiently small for its higher powers negligible 
comparison with h?, approximately given by. 


and represent the number subintervals which characteristic 

dimension the domain divided, corresponding the mesh sizes 

values and e’;, e’; are the corresponding errors, Eq. gives: 


which the extrapolated value Table gives the coefficients 


a;, the h?-extrapolation formula for the most commonly used values 
and 


TABLE 


3.769234 2.769234 


are only functions the ratios for instance, the (2/1)-extrapolation 
formula holds also for (4/2). 

If, Eq. 40a, the first two terms the series are taken into consideration 
and the formula applied three approximations k;, k;, and obtained 
subdividing the characteristic length into parts, elimination 


2 
aj 
3/2 2.28571 1.77778 
4/3 16/7 16/9 
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Table gives the coefficients for the most commonly 
common factor without changing the formulas. 


TABLE 


2.025 128/120 1.06667 
4.34028 4,096 /1,008 4.06349 
11,664/1,980 5.89091 6.31313 
5,000 /3,072 1.62760 /3,072 0.63281 
2.50104 360 1.62760 1,944/15,360 0.12656 


The formula and the formula are used 
connection with first order differences. 

Making the same type simplifying assumptions the formula for the 
error due second order differences, 


the 


Table and Table give the h‘-extrapolation coefficients and the 
polation coefficients most often needed. All these formulas are homogeneous 
and hold when all the values are multiplied the same number. 

The following points should kept mind using the extrapolation 
formulas Section 


(a) The h?-extrapolation formula and the formula assume 
that the successive approximations approach from the same side. When 
the approach the oscillating type, the formulas may used separately 
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connection with the approximations from above from below. The two 
extrapolated values thus obtained may usefully compared. unsafe 
use extrapolation without having least three approximate values 

Even the approximate values approach the true value mono- 
tonically (either from above from below), not known offhand whether 
the extrapolated value above below the true value. 


TABLE 


625 /369 256 /369 


TABLE 


3/2/1 2,187/ 1,680 1.301786 512/ 1,680 0.304762 1,680 0.002976 
5/4/3 2.170139 65,536 /50,400 1.300317 6,561 /50,400 0.130179 


(c) true that, general, the gives better results 
h?-extrapolation, but impossible determine offhand whether 
h?-extrapolation based larger values will give better results than 
the largest which the same the largest value the h?-extra- 
polation. The same statement holds true for the and the 

10. APPLICATIONS 


this section the method finite differences will used compute the 
lowest buckling load beams, plates, and shells under various load and 
boundary conditions. The following applications have been chosen: (1) Be- 
cause they illustrate the technical details and the properties the method; 
(2) because they indicate the asymptotic approach the successive approxima- 
tions the true value the buckling load; (3) because, the best 
edge the writer, the solution the problem had not been previously at- 
tempted. 

simplify and shorten the treatment, the differential equations and 
boundary conditions each problem, well the true the best available 
approximate value the buckling load, have been taken from the fundamental 
wise stated. 


(1) (2) (3) (4) 
16/15 1.06667 1/15 0.06667 
3/2 1.24615 16/65 0.24615 
4/3 0.46286 
5/4 0.69377 
(1) (2) (3) (4) (6) (7) 
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EULER’s PROBLEM 
simply supported, slender beam, length and constant moment 
inertia buckles under the action two equal longitudinal loads (Fig. 11(a)). 
Calling Young’s modulus and f(z) the normal deflection, and choosing 
the origin the longitudinal the left support, the differential equation 
and boundary conditions become!® 


and 


Fie. 


(1) First Order the length into equal parts 
length L/n and letting the ordinates the deflection curve 

Substitution the latter expression Eq. 44a gives the statement the 

boundary problem terms differences: 


Elastic Timoshenko, McGraw-Hill Book Co., Inc., New York, Y.. 
Ed., 1936, paragraph 12. 


and 


1462 
sul 
(44a) 
which the primes indicate differentiation with respect 
(a) 
- == bs 
x 
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which 


and the approximate value obtained the approximation with 
subintervals. 


Approximation 


—the first approximation and the corresponding error become 


Approximation 3.—Noticing that the deflection curve symmetrical 
with respect the middle point the beam (Fig. 11(c)) and applying Eq. 45a 


and The determinantal equation is: 


tion, which makes smaller: 1.62 and 9.54918 3.2%). 
determinantal 

and 9.64616 2.3%). 


large number successive approximations has been computed the 
foregoing example investigate the mode approach toward but 
should noticed that, since the approximations 3,5,7 involve deter- 
minants the same order the approximations 2,4,6, practice only 
the odd approximations would computed. 


Pr. 
the 
7 
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The formula (Eq. 41) with Table and the pola 
polation formula (Eq. 42) with Table give the results Table 
The extrapolated values are near the true value that practice two 
approximations and would give result more than suffi- 
ciently accurate. 
TABLE FOR THE CASE SIMPLY SUPPORTED BEAM 
witH MoMENT INERTIA 
4,5 9.86314 —0.065 4,5,6 9.86932 —0.0028 
9.86974 
(2) Second Order y’’ Eq. the boundary 
value problem (Eqs. 44) terms second order differences, becomes 
and 
which 
Since the difference equation (Eq. 48a) involves the values two points 
the left and two points the right the values and must 
known before the equation can applied and These values are 
easily determined using first order differences expressing the boundary 
conditions and these points the beam simply supported and, 
These relations state the geometrical fact that and are inflection points 
the y-curve. 
Proceeding was done before with first order differences, the following 
results are obtained: 
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the formulas and the 
polation formulas (Eqs. and Tables and 4), the following extrapolated 
values are obtained: 9.85(— 9.87(— 0.02%); and 
9.8695(— 0.0009%). 

The evaluation the second order approximations and extrapolations does 
not require increased amount labor, but the second order approximations 
are much the true value than are the first order approxima- 


tions. problems this type may well pay use second order differences 
directly. 


the previous problem, the moment inertia varies according the 
law (in which constant), 


Eq. 44a becomes, Eq. 15b, 


which 


approximations are obtained before, remembering that symmetrical 


about the middle point the beam. Table gives the results the computa- 
tions. 


TABLE FOR THE CASE SIMPLY BEAM 
WITH VARIABLE INERTIA 


h?-EXTRAPOLATIONS 


n kn en(%) 


e(%) 


16.0000 2,4 16.0000 16.5114 —0.13 
15.0000 —9.27 3,5 16.3933 3,5,7 16.4964 —0.22 
15.8917 5,7 16.4774 

16.2526 

16.1786 —2.15 


The approach one sided, but not monotonic, whereas the 
separate sequences even and odd are both monotonic. The extrapolations 
Table are obtained these sequences. 


The energy method applied this problem with assumed sinusoidal 


deflection gives upper bound 16.5342. Hence probable that the 


average the two more than 0.2% off. 


M 
I(x) 
(52) 
be 
are 
50) 
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Narrow Cross SECTION 


beam narrow rectangular cross section simply supported 
L/2, while the rotation its ends around the axis the beam pre- 
vented. The beam buckles laterally under the action vertical load 
the smaller the two flexural rigidities the beam its 
principal planes, the torsional rigidity, and the rotation section 
the boundary problem 


and 


which 


Table gives the results the computations: 


TABLE FOR THE LATERAL BUCKLING 


h?-EXTRAPOLATIONS 


11.3137 —33.2 2,4 16.21 16.905 —0.09 
18.0000 6.32 3,5 16.62 3,5,7 16.933 +0.02 
14.9827 —11.82 4,6 16.83 


The true value 16.93. this particular problem the approach 
the oscillatory type, the even being smaller, the odd being 
larger than and the extrapolations have been obtained the even and odd 
sequences. 


Elastic Stability,” Timoshenko, McGraw-Hill Book Co., Inc., New York, Y., 
Ed., 1936, 250. 
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transverse load the corresponding boundary value problem be- 


and 


and, Eq. the finite difference problem reduces 


and 
which 


Table gives the results the computations. 


TABLE FOR THE CASE LATERAL BUCKLING 
Narrow CANTILEVER BEAM 


2,3,4 4.030 +0.42 
4.014 +0.025 


The correct value 4.013 and the convergence monotonic 
tion may applied for 3,4 and 4,5, whereas the 
may applied any sequence, and particular 2,3,4 and n=3,4,5. 


SIMPLY SUPPORTED SQUARE PLATE 
UNIFORM COMPRESSION 


square plate with sides and flexural rigidity simply supported along 
its boundary, buckles under uniform compression per unit boundary 
length. Taking the parallel the plate sides with origin one 
its corners, and calling w(z,y) the deflection the plate, the corresponding 
boundary value problem becomes” 


Elastic Timoshenko, McGraw-Hill Book Co., Inc., New York, Y., 
Ist Ed., 1936, 245. 


333. 


= P?, LA _ (58) 
3.933 4,5 4.006 
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and, the boundary, 


and, the boundary, 


First Order use first order differences (Eq. and 
Fig. 4(c)) the corresponding difference problem becomes 


and, the boundary, 


which the subscripts indicate the grid points according Fig. 12, and 


Approximation 2.—Applying Eq. 6la L/2: and 
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0 4 C4 


The true value 19.739. 


SQUARE PLATE WITH CLAMPED UNDER 
COMPRESSION 


the square plate Example has clamped edges, the boundary value 
problem becomes” 


and, the boundary, 


being the normal the boundary. the first order difference operators 
Figs. 4(b) and 4(d), the corresponding difference problem becomes 


and 


which 


Table gives the results the computations. 


TABLE FOR THE BUCKLING 
WITH CLAMPED EDGES 


SQUARE PLATE 


494 


The literature the buckling rectangular built-in plates extensive, 
but rigorous solution the problem has been given. For the case the 
present example, Weinstein‘ gives very close upper and lower bounds whose 
average, 52.385, has been assumed true value for the evaluation the 


Elastic Timoshenko, McGraw-Hill Book Co., Inc., New York, Y., 
lst Ed., 1936, pp. 305, 298, and 364. 


43.342 4,5 1.58 4,5,6 53.243 +1.64 
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errors the extrapolations. (The average the bounds given Mr. Timo- 
52.455 and differs 0.13% from the average Mr. Weinstein’s 
closer bounds.) interesting notice that all the h?-extrapolations are 
from below, whereas some overshoot the mark. The 
computation the higher approximations very cumbersome, but the 
extrapolation gives excellent results the first three approximations. 


RECTANGULAR PLATE WITH CLAMPED EDGES 
UNIFORM 


The problem Example can easily extended the case rec- 
tangular plate with sides lengths and parallel and respectively. 


Eq. for 1/2 and its square, the following equation for the difference 
problem obtained: 


with 
and 
which 
TABLE FOR THE BUCKLING RECTANGULAR 
PLATE WITH CLAMPED EDGES 


137.400 153.906 

149.779 154.768 

152.972 
135.240 


The results the computations are presented Table 10.% The errors have 
been computed true value the buckling load obtained the energy 


—in which the ratio the side parallel the side parallel The 
expression for this upper bound the load 


Elastic Stability,” Timoshenko, McGraw-Hill Book Co., Inc., New York, 
Ed., 1936, 365. 


Loads for Built-In Rectangular Plates Finite Differences,” Gold, thesis submitted 
New York, Y., 1947, partial fulfilment the requirements for the degree 
aster 
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COMPRESSION 


simply supported 45° skew plate, with sides parallel and 45° sides 


buckles under uniform compression per unit boundary length. 
The corresponding boundary value problem identical with the problem 
Example 

45° skew plate may covered both rectangular and skew grid. 
Using skew grid withr and 45°, Eq. 23, Fig. Section yields 
the difference problem: 

and, the boundary, 


which 


n n2 


Table gives the results the computations. The successive values 


TABLE FOR THE CASE THE BUCKLING 


Skew CoorDINATES (b) 


polations polations polations polations 


obtained square grid are given Table together with the h?- 
extrapolations and the 

rigorous solution this problem known and value may 
accepted plausible. 

must noticed that, obtain the approximation Cartesian 
coordinates, necessary solve determinantal equation the eighth 
order. Skew coordinates are much better adapted the problem and lead 
the final result with considerably smaller amount Jabor. 


SIMPLY SUPPORTED, EQUILATERAL, TRIANGULAR PLATE 
UNIFORM CoMPRESSION 


simply supported, equilateral, triangular plate, sides buckles under 
uniform compression per unit boundary length. The boundary value 
problem solved identical with the problem Example but this 
case the plate may covered triangular grid. Eq. Fig. 9(a), 


1471 
kn kn \ 
2,3 28.73 2,3,4 27.51 16.00 2,3 25.76 2,3,4 26.68 
26.63 27.82 3,4,5 27.21 21.42 3,4 26.45 3,4,5 26.72 
27.15 4,5 27.43 23.62 4,5 26.62 
ee 
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the difference problem becomes 


and, the boundary, 


which 


and 2z, the sum the values the vertices the hexagon surrounding 
the ith point. 


TABLE FOR THE CASE THE BUCKLING SIMPLY 
EQUILATERAL, TRIANGULAR PLATE 
COMPRESSION 


(%) 
—2.64 
—1.04 +0.06 
—0.44 +0.002 
6.5 


Table gives the results the computations. The correct value 
Making use finite difference operators irregular triangular co- 
triangular plate any shape may similarly solved. 


BUCKLING SIMPLY SUPPORTED HEXAGONAL PLATE UNDER 
UNIFORM COMPRESSION 


The buckling load regular hexagonal plate, simply supported and 
subject uniform compression, may computed Eqs. and 
Example 


TABLE FOR THE CASE THE BUCKLING SIMPLY 
HEXAGONAL PLATE UNDER UNIFORM COMPRESSION 


The values the successive approximations are given Table 13, which 
the number body points considered and hence represents the order 
the determinant solved; 7.20 may taken plausible value for 


Elastic Timoshenko, McGraw-Hill Book Co., Inc., New York, Y., 
lst Ed., 1936, 371. 
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BUCKLING INFINITELY LoNG CYLINDRICAL SHELL 
AXIAL COMPRESSION 


infinitely long, thin, cylindrical shell radius thickness flexural 
rigidity and Young’s modulus buckles under uniform axial compression 
per unit length. Assuming that the shell buckles into half waves length 
the corresponding boundary value problem 


which the radial displacement. 
15, the corresponding difference problem becomes 


and 


which 


(75) 


Since the shell infinitely long, will buckle into infinite number 
half waves, whose length such make minimum. Letting 


6.066640 0806460060 4660606 04 


This value rigorously correct and obtained all the successive 
approximations the value min- 


Stability,” Timoshenko, McGraw-Hill Book Co., Inc. New York, Y., 


and 
Thus 
5 
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Eqs. and 77a, the half-wave length corresponding min is: 


the corresponding value is: 2.8284 and 9.96%. 

0.055%); 3.1370 0.15%); and 3.1468 0.17%). The 


11. CONCLUSIONS 


this paper (which essentially experimental character) the method 
finite differences applied mostly problems with known 
check the rapidity and the type convergence, well the validity the 
suggested extrapolation procedures. However, the scope the method 
such allow the attack large variety buckling problems. Among 
the advantages the method are the following: 


(a) Problems with different types boundary conditions can solved 
the same basic procedures with comparable amount labor all cases. 

(b) The variability the moment inertia beams and plates does not 
require special procedures and general does not increase materially the 
amount labor involved. 

(c) Beams with compressive loads variable along the axis the beam 
and plates with known variable distribution stress the plane the plate 
may considered without increased difficulty. 

(d) Anisotropic plates fall under the scope the method. 

(e) Plates irregular shapes may dealt with suitable set coordi- 
nates, using the equations that best fit the problem each case. 

(f) The method does not require knowledge higher mathematics. 
interesting notice that, with the exception relaxation procedures, all 
other methods solution, both approximate and rigorous, require instead the 
use integral calculus, differential equations, series expansion. 

(g) The method applicable all characteristic value problems involving 
linear differential equations, ordinary partial, with constant variable 
coefficients. 


There are also limitations and difficulties inherent the method: 


(1) The method not applicable when the distribution longitudinal 
stress the beam plane stress the plate unknown. This limitation 
common all other methods, except the energy method. 

(2) There way knowing whether the extrapolated values are greater 
smaller than the true value. 
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(3) The literal evaluation the determinantal equation for the higher 
approximations may cumbersome. 

This difficulty may avoided explained Section but interesting 
notice that all series solutions, including those used connection with the 
energy method, lead determinantal equations and hence involve the same 
amount labor the final step the computations, although sometimes 
larger amount labor required the derivation the determinant itself. 
Moreover, most the other methods not lend themselves easily extra- 
polation procedures. (Relaxation procedures, also involving finite differences, 
seem more cumbersome than the present method, since they require the 
simultaneous evaluation the deflection function, which seldom needs 
known practice. The third book Mr. Southwell relaxation methods, 
which planned cover characteristic value problems, has not been published 
the time this writing (1949) and therefore impossible make 
final comparison between the two methods.) 


The width its scope and the simplicity its application should recom- 
mend the finite difference method the attention the engineer and the 
applied mathematician. 
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STUDY END CONNECTIONS FOR STRUTS 


ASCE 


Results tests evaluate the effect different types end connections 
upon the static and fatigue strengths struts are described this paper. 
Tubes, single angles, and double angles were tested with the usual types end 
connections employed for these members well some types not common 


use. Strain measurements were made determine the distribution stress 
the members and their connections. 


was found that the rigidity the connection, degree eccentricity, and 
the amount stress concentration caused changes section the column 
ends had marked influence the static and fatigue strengths. 

Although the tests were made aluminum alloy members, believed 
that the comparisons drawn apply members other materials also. 


INTRODUCTION 


One the difficult problems the design latticed structures concerns 
the end connections web members. lightweight structures such 
towers, booms, and relatively light roof trusses, structural angles are generally 
used, and the members are connected bolts rivets through one leg. 
generally recognized that this method attachment produces eccentricity 
loading that: (1) The deformations are larger than estimated the 
basis uniform stress the member; (2) the maximum stress somewhat 
larger than the average stress; and (3) the ultimate strength less than that 
similar member with axial The performance the structure 


influenced the performance the individual members, reasonable 
comments are invited for immediate publication; insure publication the last dis- 
cussion should submitted May 1950. 
Senior Research Engr., Aluminum Research Laboratories, New Kensington, Pa. 
Research Engr., Aluminum Research Laboratories, New Kensington, Pa. 


Tests Latticed Structural Frames,” Hartmann, Moore, and Marshall Holt, 
Technical Paper No. Aluminum Research Laboratories, New Kensington, Pa., 1938. 


Hill, Technical Paper No. Aluminum Research Laboratories, New Kensington, Pa., 1940. 
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believe that the analysis the structure will simpler and more satisfactory 
the behavior the individual members more thoroughly understood. 
addition, the performance the structure might improved the eccentricity 
loading the individual members eliminated least decreased. 

The round tube with symmetrical connections appears desirable 
type member, but the difficulties making suitable end connections have 
discouraged the use this type member riveted construction. The use 
tubular members welded construction, however, has been discussed and 

experimental investigation involving both static and repeated loadings 
was undertaken the Aluminum Research Laboratories New Kensington, 
Pa., study the relative merits various types end connections for round- 
tubular struts, single-angle struts, and double-angle struts. Although alu- 
minum alloy members were used the tests, believed that the com- 
parisons drawn apply other metals also. 


DESCRIPTION SPECIMENS 


The specimens used this investigation are listed Tables and 
Types and were in. outside diameter, with walls 0.134 in. 
thick, and were seamless, drawn round tubes aluminum alloy the H32 
stabilized quarter-hard temper. Types and were single 
}-in. structural angles extruded aluminum alloy the 
fully heat treated and naturally aged temper. (Abbreviations such 17S and 
52S denote the trade designations which the alloy commonly 
Type was composed two extruded unequal-leg angles in. 
in. separated washers. The specimens were loaded 
through gusset plate connection. 

The tubes types and had flattened ends facilitate making 
riveted bolted connection. type the ends were flattened against 
tube having wall thickness 0.138 in. and length in. which was 
telescoped into the end the member serve filler. The flattening was 
done cold with two blows hammer. The resulting thickness the 
flattened end with the filler tube left was about For the static tests, 
the specimens were attached, turned studs, directly steel loading fixture. 
type there was filler tube, but filler plate was inserted temporarily 
and the end the specimen was flattened against this plate. The paddle- 
shaped end was split along the edges, and gusset plate fitted snugly into the 
slot thus formed. 

The tubes types and had slotted ends but were not flattened. 
type the ends the tubes were fitted with segmental plugs prevent flat- 
tening against the gusset plate when the bolts were tightened. type 
the tube and gussets were joined welding both sides the gusset plates. 


Saves Steel Long Engineering News-Record, July 1943, 37. 


Use Tubular Steel Sections Quisling and Cosman (rept. prepared 
under the auspices the Literature Div. the Eng. Foundation Welding Research Committee), Supple- 
ment Journal, Am. Welding Soc., March, 1943. 


Transactions, ASCE, Vol. 102, 1937, 1279, Table 
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The welds were made two passes with oxyhydrogen flame. specimen 
D5, which was tested under repeated load, the weld was continuous around the 
end the gusset plate whereas all the other specimens type the weld 
was shorter than the overlap the tube and gusset plate. 

The specimens type Table were simply equal-leg structural angles 
with the connections through only one leg, the conventional type single- 
angle connection. 


FILLER EQUAL IN THICKNESS TO THE GUSSET 
130" WALL -10 LG TUBING PLATE INSERTED BEFORE FLATTENING, 
BEFORE FLATTENING TO SERVE AS A FULLER TUBE FLATTENED BEFORE SPLITTING 


RIVETS 


GUSSET PLATE DOES NOT 
BEAR AT BASE OF SLOT 


TYPE B 


134° WALL $2S-H32 TUBING 


€wO OF TUBE FITTED SPECIMENS DI TO D4 INCLUSIVE THE GUSSET PLATE 
WITH SEGMENTAL PLUGS DOES WOT BEAR AT BASE OF SLOT. 
SPECIMEN 0S WELDED AROUND EDGE OF PLATE 


GUSSET PLATE DOES NOT 
BEAR AT BASE OF SLOT 


ee 


Lengths, inches; see diagrams. 


Tested columns with round ends. 
ends. Fatigue tested. 


Tested columns with flat 


The specimens type were equal-leg angles with flattened ends. 
prevent cracking along the fillet during the flattening operation, the heel 
the angle was rounded that the material was nearly constant thickness 
for length about 8in. The flattening was done with hammer while the 
metal was temperature about 700° The specimens were then reheat- 
treated and straightened straightener. 

The ends the equal-leg angle specimens type were only partly flat- 
tened. This operation was done with the metal room temperature using 
special dies hydraulic press. 


2M 24 
2%" WALL $25-H32 TUBING | 
hem 
7 
] 
- 
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The specimens type were double-angle members. The long legs 
the individual angles were adjacent and were separated washers 
spaced about in. apart. The spacing the stitch rivets somewhat less 
than the spacing about in. allowed the accepted method 
The gusset plates were fitted between the two angles forming member sym- 
metrical about one axis. 


TABLE 2.—ANGLE THROUGH PLATES 


A RIVETS 


TYPE F 


Tested columns with flat 


Lengths, inches; see diagrams. Tested columns with round ends. 


ends. Fatigue tested. 


The elements the cross sections are given Table 3(a). The particular 
sections used were chosen that the least radii gyration the tubes, single 
angles, and double angles would approximately equal. The mechanical 
properties the materials are given Table Both the tensile and com- 
pressive tests the tubing were made specimens the full cross section, 
and the strains were measured with type Huggenberger tensometers (gage 
length and mechanical multiplication about 2,000). The tensile tests 
the angles were made standard sheet metal specimens 0.5 in. wide,® and 


Structural Handbook,” Aluminum Co. America, Pittsburgh, 1948, 72. 


Methods Tension Testing Metallic Materials’ (E8-46), A.S.7.M. Standards, Phila- 
delphia, Pa., 1946, Pt. 1B, 297, Fig. 
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the strains were measured with Ewing extensometer having gage length 
in. The compression tests were made specimens the full cross 
section having slenderness ratio about 10. The compressive yield strengths 
were determined from stress-strain relations obtained measuring the 
relative movement the platens the testing machine dial gages the 
corners the platens. This procedure has been found satisfactory com- 
parisons such values yield strength with values determined from precise 
stress-strain curves. The values given Table 3(b) meet the requirements 


federal specification for 17S shapes the United States Navy 
Department specification 44T32 for 52S tubing. 


TABLE CONNECTIONS FOR STRUTS 


(a) Sections MECHANICAL PROPERTIES MATERIALS 


ness in 


(1) 


Round tube... 
Angles— 


0.136 0.891 0.750 24,200 17.3 


Unequal leg. 


resentative samples and the nominal specific gravity the 


Computed from the length and weight 
trade designation which the alloy commonly recog- 


material. Abbreviations Col. denote 


nized. offset. Outside diameter and wall thickness, respectively. About the axis least 


stiffness. About the axis parallel the long legs; angles are separated in. 


The nominal composition alloy 52S includes 2.5% magnesium, 0.25% 
chromium, and the remainder aluminum. The various tempers this material 
are obtained introducing certain amounts cold work after annealing 
treatment and then stabilizing aging for hours 320° for 1.5 hours 
350° The nominal composition alloy 17S includes copper and 


each manganese and magnesium. The stronger tempers this alloy 


are obtained solution heat treatment consisting heating 945° for 
short period depending the thickness the material and immediatley 
quenching cold water. The material ages naturally room temperature 
that the full temper reached about days. 


TESTING 


Round End each type was tested under static load 
column with round ends. The steel end fixtures bore against testing 
machine platens that were supported knife-edge bearings shown Fig. 
This loading condition corresponds the attachment web member 
chord zero torsional resistance. The tubing specimens and the specimens 
type were centered the platens that the axis the specimen was the 
plane the knife-edge bearings the platens. Thus, the shafts the mem- 
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bers were loaded axially although the end connection for specimen (specimen 
type was subjected some bending. The single-angle specimens were 
placed that the plane shear between the angle and the gusset plate was 
the plane the knife-edges. 

Deflections the center were measured with scale and taut wire. 
some cases the deflections points near the gussets were also measured. 


COLUMN Loap CoLUMN SPECIMEN AND ADAPTERS 
PEATED LOAD 


The static load tests were made Amsler testing machine 
maximum capacity and intermediate load ranges. The load range 30,000 
was used these tests. machine (shown Figs. and the four- 
column type and, being relatively rigid, well adapted for column testing. 
Two specimens round tubing (specimen with round ends and 
specimen with flat ends) are shown under load. The periodic calibrations 
indicate that the loads are correct within 0.5% for the load range used. 

Flat End specimens each type but different lengths were 
tested under static load columns with flat ends. The steel end fixtures bore 
against testing machine platens that were fixed against tipping and rotation 
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(see Fig. condition corresponded the attachment web member 
chord great torsional stiffness. Deflections the center were measured 
with scale and taut wire, and strains were measured number points 
type Huggenberger tensometers (mechanical multiplication about 
1,200). points rapidly changing stress, gage length 0.5 in. was used; 
other places gage length in. was used. 


Fatigue consisting pair tubes angles in. long were 
tested under repeated load ranging from zero maximum compression. 
The two-part specimen was attached the adapter the testing machine 
shown Fig. Using plate fulcra within the adapters the total load was 
divided equally between the two members. After the load range has 
adjusted, the machine was operated hand while strains were measured 
the center and near one end, again using the type Huggenberger tensometers. 
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The relative deflection the two members was measured with calipers and 
scale. 

The fatigue tests were made units and the structural fatigue testing 
the Aluminum Research Laboratories shown Fig. Essenti- 
ally each unit consists two parallel beams, one the foundation beam and the 
other the loading beam. They are hinged together plate fulcra one end 
and adjustable-throw motor-driven crank the other. The specimen 
inserted between the two beams near the hinge end and undergoes cyclic 
loading the crank rotates. The loading beam used the load weighing 
device and the deflection the beam measured while the crank turned 
slowly hand. The load the specimen then obtained from the calibra- 
tion curve previously determined. The inertia effects for the machine operat- 
ing normal speed (210 cycles per min) have been determined that the 
necessary small correction (approximately 2%) can applied the loads 
obtained from the measured deflections the loading beam. 

The machine equipped with cycle counter and each unit has cutoff 
switch. The machine stopped automatically when specimen fails. The 
arrangement sensitive that change load less than 600 will stop 
the machine. The construction the machines such that tests can made 
readily specimens only one length. 


Discussion RESULTS 


Column Strengths Static load deflection relations shown 
Fig. are typical those obtained the tests with static loads. the tests 
with round ends, the deflections points near the ends the gusset plates 
were almost great, general, the deflection the center, that the 
tubes and angles remained almost straight and the failures were double- 
hinge type (see Table 4). The specimens appeared articulated with 
joints the knife-edge bearings and the gusset plate flattening the shaft. 
should expected for members with this characteristic, the ultimate loads 
were relatively low. comparison the test results with the Euler column 
curve interest, but should borne mind that with such discontinuities 
along the length and eccentricities loading these specimens not fulfil the 
conditions Euler column. The comparison shown Fig. where the 
data points are plotted slenderness ratios corresponding length equal 
the distance between knife-edges and the radius gyration about the axis 
parallel the gusset plate the case the symmetrical sections the axis 
least stiffness the case the single angle. From this comparison, 
appears that the type specimen was most nearly like Euler column since 
supported 76% the computed load. The connection the end fixture was 
probably much stiffer this specimen than any the others since was 
attached directly the steel end fixtures. Specimen which had the 
greatest eccentricity the angle was supporting only 34% the Euler load 
(gross area) when the test was stopped because the large deflections. The 
test specimen was also stopped for this reason, but this case the load 
was 51% the Euler load. Although specimen has symmetrical section, 


Machines for Testing Structural Units,” Templin, Proceedings, A.S.T.M., Vol. 39, 
1939, pp. 711-722. 
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4in. from 
gusset plate 


LOAD, P in Thousand Pounds 


Type Type D2; 
A=68in. and A=68in. and 


Type D3; 
in. and 


MEASURED DEFLECTION, INCHES 


PLATES 


Description Over-all Maximum 
Specimen length load, Type failure and test condition 
specimen (in.) (Ib) 
(1) (2) (3) (4) 
attened end tube 
Flattened end tub 11,600 Lateral 
Tube split but ouble hinge ends.¢ 
12,250 Lateral deflection.¢ 
not flattened 18,600 Lateral 
11,900 Lateral deflection, hinged gussets.¢ 
3,500¢ Double hinges ends.¢ 
16,050 Lateral deflection and twisting.¢ 
Plain angle 16,400 Rivets 
18,400 Lateral deflection and 
29,425 Lateral deflection and 
1,500 Double hinge ends.¢ 
Flattened angle 20,150 Lateral defiection.¢ 
ouble hinge ends.¢ 
443 24,850 Lateral deflection.¢ 
2.000 Double hinge.¢ 
Double angle 84} 20,800 Gusset plate 


with round ends, with knife-edge fixtures. Tested column with fiat ends. Test stopped because 
large deflections. 


Between knife-edges. Attached directly (by turned studs) steel end fixture. Tested column 
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supported only 15% the Euler load, probably because the relatively thin 
gusset plate made effective hinge. specimen F1, the secondary enc 
hinge was the ends the flattened part the angle where the material was the 
only in. thick, compared the gusset plate. comparison 
the strengths specimens and indicates the undesirability com- 


SPECIMEN NOS. 
SPECIMEN NO.! 2 ANDO 3 


TUBING SPECIMENS 


COLUMN STRENGTH LB PER SQ IN 


3% x 3% ANGLES fs 
(F) 
SPECIMEN SPECIMEN NOS. 
2,3,4 AND 
(RIVETS SHEARED) 
ANGLE SPECIMENS 


SLENOERNESS RATIO, 


COLUMN STRENGTH LB PER $Q IN 


pletely flattening the angle. other words, desirable maintain some 

the stiffness the member the region the connection. 
Figs. and show typical failures the longer specimens tested col- the 

umns with flat ends and Fig. shows the failure the short specimens 

types and sim 


28000 
METHODS LOADING 
' 
BOLTED DIRECTLY END FIXTURES 
32000 
20000 
16000 
12000 
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Fig. will noticed that the tube appears hinged near the 
end the flattened part. The round end Euler column load for tube with 
the length equal the distance between the inside studs (70.5 in.) 10,300 
The test result 11,600 lb—showing some restraint, should expected. 


BR 


The indicated fixity coefficient about 0.94. With less rigid end fixtures 
the degree restraint, course, would have been less. 


Specimen shown Fig. 7(a) exhibits double-hinge type failure 
similar that exhibited the specimens tested with the knife-edge fixtures. 
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course, this case both hinges developed within the specimen, one the 
end the flattening and the gusset plate. Specimen (see Fig. 
with thicker gusset plate in. compared in.) did not develop double- 
hinge effect although the same two points weakness were present. Fig. 7(b) 
shows this specimen under load near the maximum. The bending the 
gusset plates apparent. appears from this comparison that the use 
rather thick gusset plates advisable with members this type, the sup- 
porting structure can furnish some restraint. 

Fig. 7(c) shows specimen after testing. obvious that thicker 
gusset plate would have increased the ultimate load this specimen. The 
computed Euler load for this member, assuming the effective length equal 
the length between inside studs the end fixtures, 25,900 lb. The test 
result 80% this value. The failure specimen H3, which the angles 
were attached directly the steel end fixtures (Table 2), was lateral bending 
parallel the plane the steel gusset plates. The ultimate load 87% the 
Euler load for specimen having effective length equal the distance be- 
tween the inside 

The plain-angle specimen (E2, Table failed lateral bending accom- 
panied twisting. Specimen (which was about in. longer than specimen 
and was attached directly the steel end fixture) supported ultimate 
load approximately 15% greater than that specimen E2. 

the shorter specimens tested with flat ends all but specimens B3, E3, 
and failed lateral deflection. The shapes specimens and after 
failure gave the appearance the double hinge The ultimate load 
specimen (which had tube length in.) was 10,930 compared 
9,550 for specimen which had tube length in. (see Table 1). 
the case the flattened angles (specimen compared specimen F2), 
the ultimate load the shorter specimen was only 85% that the longer 
specimen. Specimen failed shearing the rivets load 16,400 
The single shear stress the rivets was 26,800 per in. which 
line with the expected value for such rivets. The strength specimen 
which had the same over-all length but gusset plate was 29,425 lb. This 
specimen after failure shown Fig. The other short specimens tested 
with flat ends failed lateral deflection. the case specimen 
apparent hinge developed the flattened part, and specimens C3, D3, and 
the hinge developed the gusset plates. 

The column strengths developed the tests with flat ends are also shown 
Fig. The effective length was taken the distance between the inside 
studs the end fixtures. addition the Euler column curve, Fig. shows 
curves column strength computed the reduced-modulus and straight- 
line formulas for the tubing and angles, respectively. The reduced-modulus 
formula the same the Euler formula, except that this comparison the 
tangent modulus used the reduced modulus for values average stress 
above the proportional limit. The straight-line formula for the 
tubing, based the yield strength given Table and determined 
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and Marshall Holt, Technical Paper No. Aluminum Research Laboratories, New Kensington, Pa., 1938. 
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which P/A the column strength, pounds per square inch; the length, 
inches (between knife-edges for round end tests between inside studs 
end fixtures for flat end tests); the minimum radius gyration the cross 
section, inches; and the coefficient describing the end conditions, taken 
unity for the aforementioned lengths. The formula for the alloy," 


and the formula for the 2-in. angles" 


The fact that most the test values fall well below the curves reflects the 
uncertainty the effective lengths chosen and the effects eccentricities 
loading. The effect the initial eccentricity was magnified the cases the 
double-hinge type failure the lack rigidity the gusset plates and the 
flattened ends the members. The increase strength due rigid gusset 
plates proved the consistently higher strengths those members which 
were fastened directly the end fixtures. 


TABLE STRENGTHS TUBE AND ANGLE 
Static 


RELATIVE AVERAGE STRESS 
(Gross 


Form Description end gusset Round 
Type Flat Ends 
in. 


Tube Flattened, not slotted, bolted None 100 100 100 
Tube Not flattened, slotted, bolted 106 
Single angle olt one leg one 
Single angle Completely flattened, riveted 
Single Partly flattened, riveted 
Double angle Symmetrical about one axis, riveted 
Double angle Symmetrical about one axis, bolted None 106 


Values are approximate. Rivets sheared. 


Fig. also shows that, except for the tubular types and the strengths 
the short flat-end struts were considerably less proportion the strengths 
predicted the formulas than were the strengths the longer specimens. 
This fact can also explained partly the effect eccentric loading. The 
initial eccentricity given type connection constant, course, re- 
gardless the length the member. the member then fails the forma- 
tion the double hinges the ends with the central part the member 
remaining comparatively straight (as the case specimens and the 
strength essentially independent the length. 


nope nr 
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stress failure. 


and attached directly the end fixtures. 


show the most consistently high strength. 


(Stresses Are Compressive Unless Marked All 
Stresses Are Pounds per Square Inch) 


AVERAGE 
In.) 


Gross 
length area 


1.714 9,000} 4,985 |+5,860 


mens were 
double-angle struts. 
effective hinge. 


test. 


Center Line 


RANGE AT: 


Lateral 


further comparison the test results can found Table which the 
relative static column strengths are listed. The values are based average 
The strengths the two, type plain angle, specimens 
which were attached directly the end fixtures, are seen 84% and 80% 
the strengths the corresponding type tubes flattened over filler tubes 


the specimens attached the 
end fixtures through gusset plates, the type tubes with segmental filler plugs 


Ratios 


11,640 
14,160 1,440 


15,200 5,540 
8,820 6,210 
9,980 4,730 
9,230 4,560 

11,490 5,240 
12,740 6,230 
9,780 1,990 
10,470 4,170 
8,300 +13,950 
8,190 9,390 
§,360 +10,360 
7,650 +11,060 
7,760 4,260 


Load which strain measurements were made, 


tions when compared with the single-angle struts. 


Measured Stresses, Static measured strains were interpreted 
into stresses simply multiplying the modulus elasticity; thus, for 
10,200,000 per in.; and, for 178, 10,600,000 per in. 
The measured stresses thus obtained are summarized Table and specific 
examples are shown Fig. The dotted lines indicate the value the 


Elastic Constants for Wrought Aluminum Alloys,” Templin and Hartmann, 
Technical Report No. 966, National Adviscry Committee for Aeronautics, Washington, January, 


should remembered that the gusset plates for the single-angle speci- 
in. thick, whereas gussets were used for the tubular and 
Obviously the thinner gusset plate should the more 
Weakness the gusset plate probably accounts for the low 
comparative strength the type double-angle specimen the round end 
the other hand, the flat end tests double angles proved 
superior single angles, despite the failure the gusset plate specimen H2. 

The data Table indicate generally greater structural efficiency for 
the tubular and double-angle types struts with relatively rigid end connec- 


1 
Spec Loads End tion 
(Ib) Mea- center 
pute Maxi- Maxi- in. 
at; line mum i mum fe 
0.904 500 10,510 10,520 9,380 18,460 0.034 
20,690 0.294 
28,380] .... 
0.130 
9,970 0.024 
12,330 0.40 
9,870 0.260 
11,860 0.090 
17,220 0.032 
12,950 0.021 
18,040 0.042 
9,540 0.030 
in. 
ths 
ths 
re- 
ber 
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average stress the gross section strut. Compressive stress plotted 
the right and tensile stress the left. The loads which the measure- 
ments were made corresponded average stresses ranging from 5,250 per 
in. for the single angles 10, 510 per in. for the tubes types and 
Although there considerable variation measured stresses, the averages 
the stresses measured the center agree very well with the simple (P/A)- 
stresses. thirteen the seventeen specimens the average measured stress 
within the computed average stress and only the case specimen 


Scale, pounds per square inch 


the difference greater than 8%. general the range stress the 
end greater than that the center, and for the specimens with flattened 
ends much greater. The maximum compressive stresses are also generally 
greater the end than the center. For example, the case specimen 
under load 9,500 the stresses the center range from 9,380 per 
in. 11,640 per in. compression whereas the end the stresses 
range from 2,850 per in. tension 18,460 per in. compression. 
The maximum compressive stress the end 1.6 times the maximum com- 
pressive stress the center. Because large lateral deflections, the range 
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measured stresses the center was greater than the end for specimens 
C2, E2, and E4. 

study Fig. indicates that, general, the stress distribution the 
flattened ends the tubes alike, despite the difference the attachment 
the gusset plate. can noted that tensile stress was developed the 
tapered part the tube near the flattening. The stress gradients this 
region are rather high. For specimens types and the stress gradients 
are relatively flat. 

the plain-angle specimens (type the load eccentrically applied, but 
addition the eccentricity there are certain restraining moments that de- 
crease the effects the eccentricity. The restraining moment one direction 
rises virtue the bending stiffness the gusset plate and the other 
direction the use more than one rivet the connection. working 
backward from the measured stresses, the actual eccentricities the flat end col- 
umns were found be, approximately: normal the gusset plate 0.75 in.; 
and parallel the gusset plate 0.31 in. 0.60 in. The nominal eccen- 
tricity taken equal the distance from the back the angle the centroid 
the section 0.94 in. The stiffness the gusset plate accounts for the 
difference. 

The eccentricities single-angle connections were greatly reduced 
the flattened ends types and The measured stresses indicated 
effective eccentricity normal the gusset plate approximately 0.09 in. 
toward the heel the angle type and 0.16 toward the toes the angle 
type with the ends only partly flattened. Thus, stresses measured the 
center the specimens types and were all compressive stresses, whereas 
the center specimen the stresses varied from 11,490 per in. 
compression 9,800 per tension. appears that greater amount 
flattening the G-type specimens would have reduced the eccentricity. 

the case the flattened end tubes, deforming the ends the angles 
produced high stress gradients near the connections. Both maximum com- 
pressive and maximum tensile stresses measured the ends the specimens 
types and are higher than those found the specimens type 

Placing two angles back back make the strut symmetrical about 
one axis greatly reduced the stress variation. The range stress the 
center specimen was only from 6,080 per in. 7,760 per in. 
compression. 

Results Fatigue distributions the fatigue specimens were 
found quite similar the distributions found the static tests flat 
end specimens shown Fig. 

the case the angles, two specimens each type were tested fatigue, 
one with the same average stress that imposed the tubing (approximately 
7,500 per in.) and the other with the same total load. the second 
specimen each type the stress cycle was from zero approximately 4,000 
per in. 

Table summarizes the results the fatigue tests giving ratios mea- 
sured stresses and the number cycles compressive loading imposed 
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each specimen. Three the specimens—C4, and H5—withstood approxi- 
mately 25,000,000 cycles stress without visible signs failure. For the 
other specimens, the number cycles produce visible cracks ranged low 
71,900. this case (specimen E6), failure occurred shearing the rivets 
which the average shear stress was about 20,400 per in. The failure 
specimen also involved the failure rivet. The average shear stress 
the rivet was about 15,300 per in., and the number cycles failure 


Stupy Enp FoR STRUTS 
(Stresses Are Compressive Unless Marked 


Srress per In.) 


Average Maximum 
Remarks 
Measured 


Measured 
P/A Center Line 


(9) 


13,650 8,000 12,700 


13,520 7,710 9,070 Removed 549,800. 
13,700 8,110 9,180 failure. 
Small cracks ends gusset 


plates; 25,208,800 cycles 
13,600} 7,5 7,610 8,650 8,690 |22,891,100 extended over 


fourths the ‘circumference. 


13,760 7,690 8,390 9,130¢ ends gusset 


26,370 7,550 Rivets sheared. 


25,610 443,800 

13,640 3,860| 5,508,200 


26,600 6,910 11,300 4,216,800 
13,580 3,820 6,080 4,950 failure. 


Maximum load stress cycle and load which strain measurements were made. Twoeach. Stresses 
ends diameter parallel the gusset plate average 9,040 per in. compared average 
7,888 per in. for the two diameters angle 45° the gusset plate. specimen D4, the aver- 
ages for all the diameters were practically the same. 


was 1,118,600. These two test results lie well above the (S-N)-curve indicated 
the data for joints with cold driven rivets reported Hart- 
mann, ASCE, Lyst, and Andrews, Jun. ASCE, This 
difference expected since the riveted joints failed tension the 
plates whereas these specimens the failures were the rivets. 

The test results Table further indicate the superiority symmetrical 
connections over unsymmetrical connections and show that flattening the ends 


Hartmann, Lyst, and Andrews, Advance Restricted Report No. 4115, National Advisory Com- 
mittee for Aeronautics, Washington, C., September, 1944. 
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tubylar members decreases their resistance repeated loads. Photographs 
some the specimens are shown Figs. and 11. Because the local 
nature the fatigue failures, expected that the strengths and lives 
the specimens were little affected the length the specimen. 

General Remarks.—The results the tests described this paper corrobo- 
rate the general conclusion Van den ASCE, based com- 
pressive tests single angles with different types connections—namely, 
that the compressive strength angles affected largely the eccentricity 


the applied loads. Although Mr. Van den Broek found that angles with 
completely flattened ends showed greater eccentricity and less column strength 
than angles connected only one leg, the opposite effect was found the 
tests described this paper. Since the eccentricity the connections 
flattened end angle depends the manner which the angle flattened and 
the face which the supporting fixtures are attached, this comparison cannot 
further analyzed because the details the end connections are not dis- 
cussed Mr. Van den Breek’s paper. 


Connections and Rivet Holes Ductility and Strength Steel Angles,” Van 
den Broek, Civil Engineering, February, 1940, 83. 
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Various types welded connections for tubular struts were described 
Quisling and 1943. With reference flattening the 
ends tubes they conclude that tubes with flattened ends 


suffer from the weakness inherent drastic changes the shapes 
the members, always objectionable. The majority them are also 
less rigid laterally, lacking the stabilizing effect the tubular shape the 
joints, where most 


These statements seem borne out the results these tests the 
Aluminum Research Laboratories. However, the comparatively high static 
strengths the type flattened end tubes (Table bolted directly the end 
fixtures indicate that the loss ridigity due flattening the ends may 
somewhat lessened first inserting filler tube reduce the degree flatten- 
ing and using very rigid gusset plates. 


REMARKS 


Although the tests made for evaluating the various types end connections 
were made aluminum alloy members, believed that the general con- 
clusions drawn are applicable members other materials. 


Round tubes are superior static column strength single angles, 
provided that the end connections are sufficiently rigid. the types 


specimens that were fastened the end fixtures through gusset plates, the 
type unflattened tubes, slotted and bolted through segmental filler plugs, 
showed the most consistently high strength (see Table 5). Based the aver- 
age stresses failure for the types and specimens which were attached 
turned studs directly the end fixtures, the static strength round tubes 
type with the ends flattened over filler tubes, was about 20% greater than 
the strength single angles the same length, fastened through one leg. 

Double-angle struts, symmetrical about one axis, are superior static 
column strength single angles, provided that the end connections are 
ciently rigid (see Table 6). The static strengths the various long single- 
angle struts tested columns with flat ends were from about 76% 96% 
great the average strength the double-angle strut tested column 
with flat ends. 

The stiffness the gusset plates and the rigidity the connections 
are important the static strength struts. 

The static strengths the round tubes tested columns with round 
ends were from 24% 76% the values computed the Euler column for- 
mula taking the effective length the length between knife-edges the sup- 
porting fixtures. The corresponding range for the angles was from 15% 

The strengths long struts the types tested, loaded columns 
with flat ends, ranged from 20% below 27% above the values computed 
the Euler formula using fixity coefficient equal unity and effective 
length equal the length between inside studs the end fixtures. 


¥ 
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The strengths the short tubular struts tested columns with flat 
ends ranged from 61% 108% the values computed the straight-line 
column formula taking the effective length equal the distance between inside 
studs the end fixtures and fixity coefficient equal unity. The corre- 
sponding range for the single-angle struts was from 37% 55%. 

the types specimens tested, double angles and round tubes bolted 
welded gusset plates placed slots the unflattened ends showed the 
most favorable distribution stress near the connections. The ratios 
maximum measured stress average stress for the various types connections 
were follows: 


Description end connection Ratio: 
Slotted tube, unflattened, double angles................. 1.38 
Flattened tube, not slotted, filler tubes ends.............. 1.75 
2.74 
Single angle, fastened through one leg..................... 3.12 


The stresses were measured loads ranging from 7,500 12,300 de- 
pending upon the type specimen (see Table 6). 

the types tubular specimens tested, round tubes with the ends 
unflattened and bolted gusset plates placed slots the ends showed the 
greatest resistance repeated loads producing stress cycle from zero about 
7,500 per in. compression. Similar specimens welded rather than 
bolted the gusset plates showed slightly less fatigue strength, whereas 
flattened-end tubes showed the lowest resistance repeated stress. 

The fatigue lives the various tubing specimens tested were follows 
(cycle average stress from zero about 7,500 per in. 


Description end connections Number 
stress 

25,004,500 

Slotted ends, bolts through segmental filler plugs............... 
22,891,100 

Slotted ends, welded gusset 20,785,400 
Flattened ends, not slotted, unsymmetrical connection............ 342,200 
Flattened ends, slotted, symmetrical connection.................. 172,000 


Angles with the ends flattened order load the member close the 
neutral axis showed greater fatigue strength than angles connected one leg, 
but less fatigue strength than double angles symmetrical about one axis. The 
fatigue lives (number cycles) the various angle specimens tested, with 
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average compressive stresses from zero about 4,000 per in. and 7,500 
per in., respectively, were: 


Description end connections per in. per in. 
Double angle, symmetrical connection. 4,216,800 
failure 
25,489,300 
Completely flattened angle.............. 1,118,600 
Partly flattened 5,508,200 443,800 
Single angle fastened through one 4,276,700 


Rivet failure 


10. desirable combination static strength and resistance repeated 
loads furnished round tubes, unflattened, with gusset plates placed 
slots the ends, and double angles, symmetrical about one axis. 

11. Flattened end tubes are roughly comparable static strength and re- 
sistance fatigue single angles fastened through one leg. 
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DISCUSSIONS 


PROPERTY SURVEYS MUST FIT THEIR TITLES 


Discussion 


WILLIAM WATTLES 


and discussers this paper should 
realize that any rational enlargement the subject matter, beyond the basic 
principles and generalities given, would produce volume far beyond any 
reasonable scope. example this scope afforded the valuable work 
this subject, presented the late Ray Skelton, ASCE, 

Each survey land boundaries establishing the location record title 
and relating the same physical conditions presents particular problem, 
different most its aspects from any other, except basic principles. 
These differences are responsible for many the variances legal decrees 
given superficially similar cases. 

Mr. Borg’s claim that certain statements the paper may misleading, 
quoting legal citations support his discussion, not proved. The paper, 
taken whole, indicates the uncertainties present any survey, and cautions 
against the flat acceptance any principle itself, but requires its analytical 
association with others. Item (‘‘Any other matters connected with the 
property the list defining (determining) boundaries, includes 
examination any citations covering matters similar the problem hand. 

The question reliability superiority physical monuments over the 
description calls record matter investigation and analysis. the 
surveyor making location from record, general rule give secondary 
value physical monuments whatever kind, against the calls and non- 
physical legal monuments recited the description. the description calls 
for physical monuments, those existing the ground can sufficiently 
identified position with the record, then the monuments will have superiority. 
subsequent court decision, adverse this rule, should construed 
modification because other factors the case, and not failure the rule. 
Mr. Borg correct stating the intent_of the parties will finally 


paper William Wattles was published February, 1948, Proceedings. Discussion 
this paper has appeared Proceedings, follows: May, 1949, Robert Borg, and Bauer; 
June, 1949, Rupert Carroll; and September, 1949, Henry Southerland, Jr. 
Glendale, Calif. 
Legal Elements Boundaries and Adjacent Properties,” Ray Skelton, Bobbs-Merrill 
Co., Indianapolis, Ind., 1930. 
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govern, but the court, not the surveyor, who may assume intent 
beyond that specifically expressed the record. 

Mr. Bauer’s discussion excellent and complete summary the neces- 
dary coordination between survey and title. Also, points out the methods 
attaining such coordination through education the surveyor and layman 
recognize facts opposed assumptions. 

Mr. Carroll complimented his discussion survey difficulties 
Texas; valuable addition survey information. His statement re- 
garding knowledge the practice previous surveyors very pertinent and 
applies elsewhere well Texas. Often such knowledge marks the dif- 
ference between correct and incorrect location. 

impossible answer fully Mr. Southerland’s implied questions within 
the scope this discussion, for approaches the realm original error 
record titles and the corrective procedure thereof; this not within the au- 
thority the surveyor nor the province the title location survey. 

perfectly true that many descriptions are not from their 
own statements recitals and that such descriptions, conveyance 
any kind, cannot pass perfect title; therefore, any survey under such condi- 
tions cannot fit title, but must considered the light fact-finding 
procedure for purposes adjustment legal methods. 

recognized that surveys are made for many purposes other than 
property boundary relocation. such cases the matter relative values and 
accuracies may enter. property survey accord with the title must 
particular and exact such title can determined; any departure, such 
assumption possession boundaries and agreed lines not confirmed record, 
converts the title location survey another form. 

The purpose the paper was present series principles property 
boundary location founded the record title the best and most consistent 
base operation. These principles must associated modified need 
be, any particular description case. the event inability 
establish correct technical position, because error insufficiency 
the record, the problem resolves into pure location for purposes adjust- 
ment legal methods. 

Matters cost, relative values land and survey expense, mechanical 
technique, allowable errors, and the like, have part the analysis property 
survey conformable with title. 

The writer wishes thank the discussers for their timely, pertinent, and 
friendly consideration the paper. 
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DISCUSSIONS 


AIR-INLET VALVES FOR STEEL PIPE LINES 


Discussion 


size air-inlet valve required maintain certain 
internal pressure pipe line during transient flow conditions can deter- 
mined easily the method proposed this paper. Undoubtedly, there are 
many conditions that may exist pipe line which these data would apply, 
and one these special cases that will described herein. 


air inlet valve 


Maximum canal 
water surface 


the discharge outlet the pumping plant pipe lines Grand Coulee 
Dam Washington there siphon shown Fig. The problem 
breaking the siphon action for return flow through the discharge line following 
power failure the pumps arose. This problem was not associated with the 
collapse the pipe line but rather with the loss water from the canal 
siphon action maintaining return flow through the pump. 

The siphon head available was limited and depended upon the canal water 
surface elevation, relative the siphon invert elevation. This relationship 
determined the internal pressure that had obtained and maintained 
the discharge line break the siphon action. 

addition this internal pressure, the rate which air must supplied 
had determined. This aspect the problem was handled recourse 


paper John Parmakian was published June, 1949, Proceedings. 
4Engr., Branch Design and Constr., Bureau Reclamation, Denver, Colo. 
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the pump characteristics the model pumps for Grand From these 
characteristics the steady-state reverse discharge normal head was obtained 
approximately 0.5 times normal discharge. Transient values 
reverse discharge were higher than this steady-state value but were considered 
too short-lived basis for selecting the size air valve needed. 

Then, knowing the internal pressure necessary prevent siphon action 
and expressing the steady-state reverse discharge terms velocity, the 
curves shown Fig. were used obtain the size the air value required. 
was concluded that air valve in. diameter, its equivalent, was 
needed. air valve would supply air rate equal the steady-state 
reverse discharge through the pumps with the internal pressure the pipe line 
the valve adequate break the siphon action. 


A.S.M.E., Vol. 63, 1941, 251 


take: 

pare 
situs 

the 

crea 

thos 

less 

sult 
sun 
slu 


AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


DISCUSSIONS 


SOME EFFECTS ANAEROBIC DIGESTION 
SEWAGE SLUDGE 


Discussion 


taken the time analyze this paper carefully, and have advanced conclusive 
thought and criticism regarding the experimental work. response these 
and the many written and verbal comments made directly the writers 
desired state that the paper was not designed advance plea for special 


method sewage sludge disposal but rather suggest tool heretofore ap- 
parently not investigated and one which might find logical application certain 
situations. 

Mr. Abbott directs attention the fact that returning digested sludge 
primary effluent for disposal the volatile suspended solids, indicated 
the writers’ computations, may increased from total 30% total 
65%. His position well taken, but fails indicate that the 35% in- 
crease volatile solids would result from the return the primary effluent 
those volatile solids the digested sludge which are slow digest and which, 
general, are nonputrescible and have low 5-day biochemical oxygen de- 
mand The volatile solids returned with the digested sludge are 
similar composition leaf mold, peat, some such matter. 

Mr. Abbott’s reluctance apply the 5-day B.O.D. test sewage sludges 
not unusual, although believed would firmer ground were 
less reluctant apply the test digested sludge than crude sludge. Re- 
sults the test applied sludges are erratic times; but, for what the 
results may worth, they are likely measure, indeed, overmeasure, the 
oxygen demand placed upon flowing stream, lake, salt water disposal 
site digested sludge which will not strand. Mr. Abbott appears as- 
suming that the remaining relatively stable volatile matter the digested 
sludge will form sludge banks. The writers’ admonition stated the paper, 

paper Rawn and Candell was published November, 1948, Proceedings. 


Discussion this paper has appeared Proceedings, follows: May, 1949, Abbott, and Robert 
Leaver; and September, 1949, Genter. 


Chf. Engr. and Gen Mgr., Los Angeles County Sanitation Dists., Los Angeles, Calif. 
Denver, Colo. 
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not return the digested sludge primary effluent for disposal there 
danger serious stranding. 

Mr. Leaver questions the desirability removing crude sludge from the 
sewage and then returning digested form for disposal. His criticism 
appears stem from the fact that digested sludge valuable natural re- 
source and should not wasted the degradation streams. the latter 
regard, Mr. Leaver directs attention the accumulation crude sludge 
eddies the rapidly moving stream waters the State Washington. Un- 
fortunately, tends apply his reasoning the Columbia River below 
Portland, Ore., into which river the crude sewage from cities tributary 
streams and the banks the Columbia has been emptied for years without 
significant effect the waters that great stream. The effective reduction 
the readily putrescible content large proportion the suspended solids 
the crude sewage Portland can surely have only beneficial effect 
existing conditions; and, applied the Portland situation, the writers believe 
that the return digested sludge the primary effluent and its discharge into 
the Columbia River will escape detection even the best laboratory 
techniques. 

Again the writers state that, the solids will strand, the practice re- 
turning digested sludge primary effluent may not satisfactory and 
probably should practiced only experimentally and with caution until its 
merit demonstrated. gratifying note that the civil engineering 
department the University Washington, Seattle, plans extend the 
writers’ research. 

The writers agree with Mr. Leaver that digested sewage sludge usually 
excellent soil builder and that nothing should allowed curtail the 
activities those who desire reclaim and utilize this material for soil building 
course, that the cost such reclamation borne those who 
profit from the venture. Surely, nothing the writers’ discussion can 
construed indicate that they (who are among the most enthusiastic pro- 
ponents the return digested sludge the soil for agricultural betterment) 
desire waste this material when can utilized properly and economically. 

Mr. Genter questions the validity the conclusions drawn from labora- 
tory experiments digested sludge. However, his position not well taken 
with respect laboratory digestion that designed illustrate the changes 
sludge allowed continue. digestion without further accretion 
crude sludge. This condition approached the later stages controlled 
stage digestion currently practiced, and demonstrated 
the latter part the curve controlled laboratory digestion. Mr. Genter 
states that: the actual plant system fresh sludge added daily the 
digesters This statement true only the primary stage; not 
applicable the later ripening stages modern digestion. 

Mr. Genter’s comments the concentration solids the crude sludge 
not apply the conditions described the paper. The solids content 
the crude sewage the controlling factor and not the concentration the 
solids the crude sludge. appears logical that the greater the strength 
and solids content the original sewage the less the primary effluent therefrom 
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will affected percentage-wise the return digested sewage sludge which, 
having originally contained nearly 40% the volatile suspended solids the 
sewage, contains less than half that amount after some days digestion 
and only 10% its originally measured 5-day B.O.D. 

Mr. Genter’s the writer believes that this highly tele- 
scoped method sewage and sludge disposal can serve only postpone the 
day real not deemed appropriate the authors. They 
feel that the recommendations contained the paper will result the final 
disposal digested sludge completely and adequately were incin- 
erated. Hundreds thousands tons similar organic material from farms, 
forests, and meadows are disposed annually inthesame manner. Again, 
response Mr. Genter’s criticism regarding the wasting this valuable 
fertilizer, Mr. Rawn desires state that the sewerage authority under his 
direction has produced and sold more digested sludge the form fertilizer 
than any authority west Chicago, While still maintaining that the 
return this organic material the soil great value, insists that there 
are locations where cannot even given away and where, processed 
soil nutrient, the operation should the expense those who will 
profit the operation. Unless the manufacture fertilizer necessary 
step the operation sludge disposal, such locations the cost should not 
borne the public. 
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DISCUSSIONS 


STREAM FLOW VARIABILITY 


Discussion 


VEN CHOW, AND DANA Woop 


Ven the computation procedure proposed the 
paper, the variability index can expressed the formula: 
= (y 9) (1) 


which the logarithm the selected discharge 10% interval the 
duration curve; the mean the deviation from the mean; and 
the number selected discharges, the given procedure. Ina 
practical operation, seldom convenient take the deviations from the 
actual mean, since such deviations usually involve decimals which are cumber- 
some handle when squared. For greatest convenience (especially machine 
computation), another form Eq. obtained utilizing the mathematical 
principle that the sum the squares deviation from the mean equal the 
sum the squares deviation from zero, less the product the total and the 
mean. Modified slightly, this principle expressed 


Eq. known the Gaussian formula, which used estimate the standard 
deviation population statistics. Eq. the mean the squares 
The variability index can defined mathematically and simply this 
formula. 

The numerical method computing the variability index was primarily 
presented for the purpose eliminating the effects introduced variations 
personal judgment. However, the method used the values discharges 
read from the duration curve plotted, did not eliminate the personal factor 
introduced plotting the curve. complete elimination the influence 
personal judgment possible only replacing the mathe- 
matical process—that is, the values and contained Eq. should 


paper Lane and Kai Lei was published September, 1949, Proceedings. 
Hydr. Eng. Laboratory, Univ. Illinois, Urbana, 
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computed directly from the raw data. Such data can grouped the 
method statistical classification reduce the number data handled 
computation. necessary, the computed standard deviation can corrected 
applying the Sheppard correction. 

Since the mathematical part the paper has been omitted, the method 
producing duration curve from the given values the variability index and 
the mean annual runoff cannot checked readily. However, the stipulation 
(under the heading, ‘‘Use the Variability Index’’) that— 


the ratio the discharge exceeded 15.87% the time the 
discharge exceeded 50% the time was equal the variability index 


—would make the variability index greater than unity for most streams, because 
the discharge that exceeded 15.87% the time usually greater than the 
discharge exceeded 50% the time. Apparently this estimate does not agree 
with the values the variability index listed Table 

The knowledge hydrologic statistics has been developed and extended 
greatly recent years. Many have found that the distribu- 
tion most hydrologic phenomena skew. Probability plottings the 
Hazen type were found unsatisfactory for skew frequencies. The computation 
asymmetric probabilities now possible the field hydrology; and the 
development special probability opened the way for straight- 
line plotting the duration curve. Therefore the study stream flow 
would advanced suitable asymmetric probability plotting could 
adapted it. The values discharge should plotted linear scale, 
that the variability index would equal the standard deviation, 
statistical sense. 

Many hydrologists like use the coefficient variation, instead 
the standard deviation, for the study variability. The coefficient 
variation equal the standard deviation divided the mean, 


For purposes comparison, the coefficient variation preferable the 
standard deviation, because the former, being dimensionless, eliminates the 
effects the difference from the mean. This elimination desirable, especially 
for stream discharge more arid regions where the annual rainfall tends 
smaller well more variable. 

For any statistical study, the quantity the available data important 
the quality. long the quality kept homogeneous, the longer record, 
larger number data, would make possible more accurate estimates the 
frequencies stream flow. the Soil Conservation Service, 


Asymmetric Probability Function,” Slade, Jr., Transactions, ASCE, Vol. 101, 1936, 
35. 


Estimated Probability Method,” Gumbel, Engineering News-Record, June 14, 
» DP. Ue. 


Line Plotting Skew Frequency Goodrich, Transactions, ASCE, Vol. 91, 
»p. i. 


Simple Method Estimating Flood Ralph Powell, Civil Engineering, Febru- 
ary, 1943, 105. 


the Gumbel Method for Computing Probability Curves,” Potter, 
No. SCS-TP-78, Soil Conservation Service, U.S.D.A., Washington, C., May, 1949. 
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United States Department Agriculture, has made probability study 
rainfall and runoff. The study was confined annual maximum values; but 
worthy note that the frequency distribution function years 
record, which measure the number data used. Since the variability 
index measure the slope the duration curve frequency curve), 
should some function the length record. Therefore, the variability 
index would more statistical value the two were compared the basis 
the same length record. 


Dana ASCE.—The selection analytical method 
depends upon the use which the results are put and the time available 
for the study. Cost also factor any engineering office. Two suggest- 
ions concerning the duration curve may value those desiring extend 
this work. 


(1) Instead using scale cubic feet per second per square 
mile,” substituting the average flow” for the period record under 
investigation will result relative variability curves flow. The values 
departures 10% intervals can then obtained similarly from this curve. 
This procedure eliminates the size the drainage area and the magnitude 
the runoff any comparisons. stream having flow 0.5 per sec 
per mile can compared with one having flow 2.5 per sec per mile 
with respect the slopes the duration curves—that is, the magnitude the 
departures from the average, regardless the size and other characteristics 
the drainage area. The writer would like know this approach might 
any way simplify the logarithmic calculations. 

(2) The second point importance that the same period years should 
used when comparing different records the same general locality. 
well known that the magnitude the runoff varies cycles wet and dry 
years, although little enough known about the amplitude and magnitude 
the curve variation. However, the Tennessee basin and numerous other 
basins, the years from 1930 1948 have been generaily the part 
cycle. The previous nineteen years, from 1910 1929, were generally 
flat part the cycle having normal annual variations. quite obvious 
that any analysis could lead erronous conclusions records the two sepa- 
rate periods were compared. This point well illustrated the doubt 
regarding the analysis records different periods years when attempting 
determine the frequency occurrence floods great 
However, for more distant localities, where the runoff cycles coincidentally 
may radically different, further clarification and the development methods 
adjustment are needed. 


This paper offers suggestions for many interesting exploratory studies, 
especially for the younger engineers training these fields who are not afraid 
undertaking many laborious calculations. 


Chf., Power Studies Branch, Div. Water Control Planning, TVA, Knoxville, Tenn. 
and Probable Future Floods,” William Creager, Civil Engineering, November, 1939, 
668. 
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DISCUSSIONS 


CONTROL THE HYDRAULIC JUMP SILLS 


Discussion 


MICHELE VIPARELLI, RALPH WEAVER, AND 
SHERMAN WOODWARD AND FLOYCE YATES 


MIcHELE this remarkable study, the authors have devised 
set simple rules for determining the height, sill that will create 
true hydraulic jump selected point rectangular channel, the flow 
being supercritical. Observations two types sill are reported—a sharp- 
crested weir and vertical jog the channel bottom. the two, the sharp- 
crested weir seems least practical because, more commonly, such barrier 
would the form broad-crested weir. Convenient graphs are pre- 
sented for the solution equations developed from laws that govern the 
creation the jump, and from laws that apply fluid flow over weirs and 
open channels. The use nondimensional coordinates makes the application 
the charts exceptionally simple. 

The same problem proportioning the sill can solved very simply, 


without recourse design graphs, with the following formula proposed 
Guido 


2 2 2 


2 
Kq. can applied directly for the case the abrupt 


rise, which the jog the channel generally created excavation the 


paper John Forster and Raymond Skrinde was published April, 1949, 
Proceedings. Discussion this paper has appeared Proceedings, follows: September, 1949, En- 
Yun Hsu, and Steponas Kolupaila; and October, 1949, Charles Jaeger, and Louis Laushey. 


Asst. the Inst. Hydraulics, Univ. Naples, Naples, Italy. 


proporzionamento traverse soglia prolungata per dissipatori Guido Nebbia, 
Acqua Gas, Roma, 1938. 


dissipatori salto Bidone, Basi Guido Nebbia, L’Energia Elettrica, Milano, 1940. 
dissipatori salto Bidone, Ricerca Guido Nebbia, ibid, Milano, 1940. 
dissipatori salto Bidone Norme proporzionamento,” Guido Nebbia, ibid, Milano 1940. 
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original bottom depth between sections and 3-3 (see Fig. 18). 
Then the designer generally knows the value the unit flow the depth 
and the depth which dependent the conditions upstream from section 
0-0. Fig. shows that: 


and 


and, then, from Eq. 39b. 


NY, 


the case weir (sharp crested broad crested) the solution for 
trial and error, solving for successively assumed values Since the 
value changes relatively little with corresponding values the final 
solution generally requires more than two trial computations. 

Kq. has advantage over Fig. which depends sharp-crested weir 
coefficients, and therefore convenient when the weir has broader other 
structurally resistant crest. 

their interpretation the experimental results, the authors properly 
state that, had the length been sufficient develop the jump completely, 
the experiments would have checked the theory more exactly. Mr. Nebbia’s 
confirm this observation. Under the heading, Ap- 
plication Results,” the authors propose that the curve for 
Fig. used, under conditions maximum Mr. Nebbia has 
presented evidence that there are conditions which the highest value 
can found discharges less than the maximum, especially when the total 
head, the flow lower stage constant under varying 


Weaver.**—In many instances has been impossible use 
hydraulic jump because either the tailwater may have been too low, which 


Hydr. Engr. Hydr. Laboratory, Norris, Tenn. 
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case the jump would not occur the desired reach, may have been too 
high, causing the jump submerged. The authors have made extensive 
study two solutions for the first condition. each, the water depth was 
raised, weir action, the depth required form jump. 

complete theoretical analysis, based the assumption that submergence 
the weir would not occur, was made for weir and abrupt rise the 
channel floor. The formula for unsuppressed weir with velocity 
approach was used with corrections applied account for suppression 
the weir, velocity approach, and nonaeration the nappe. The theoretical 
analysis interesting but, inasmuch the values the correction coefficients 
used the derivation are approximate and the inflexibility the formula 
limits its application, the practicality the equation questionable. 

The results the test agreed reasonably well with the theoretical analysis, 
but discrepancy existed which may explained tests that were con- 
Jun. ASCE, and the Carnegie Institute Technology, Pittsburgh, 
Pa. These tests were conducted determine the effect sill the tail- 
water depth the hydraulic jump. The curves Fig. indicate the presence 
velocity distribution such that which existed the foregoing tests. 
When the ratio was small, the discrepancy between the theoretical and 
the actual curves was large, whereas the theoretical and actual curves were 
almost identical when was relatively large. the larger values 
the velocity distribution became more uniform than any previous value 
The theoretical formula applied actually only the high values 

The were conducted 40-ft flume square with adjustable 
inclined weir permit fine adjustments the tailwater depth. Sills in. 
wide, varying heights extending the full width the flume, were mounted 
such manner that the force the water impinging the sill could 
measured standard electrical strain gage methods. The initial depth was 
controlled allowing the water flow under inlet plate that was fixed 
0.2 above the floor the flume. The general testing procedure consisted 
measuring the force the sill various locations throughout the jump and 
the tailwater depths the end the jump. The results from the tests 
which the end the jump was over the sill were plotted the dimensionless 
ratios the initial depth, the critical depth, and the measured 
tailwater depth, the calculated theoretical tailwater depth, required 
form jump. was found that the ratios ranged between 85% 
and 100% and that the values were function the ratio and 
the sill height. The value decreased the sill height increased. the 
sill height increased, height was reached which jump would not form 
upstream from the sill. 

was also found that the force the sill decreased rapidly minimum 
the end the jump was moved upstream position approximately over 


Rectangular End Sills the Hydraulic Jump, Level Rectangular 
Laushey, Carnegie Inst. Technology, Pittsburgh, Pa., 1947. 

Rectangular End Sills the Hydraulic Jump Level, Rectangular Channel,” 
Anderson, Carnegie Inst. Technology, Pittsburgh, Pa., 1948. 

Study Forces Rectangular End Sills Hydraulic Jump Level, Rectangular 
Hallenburg and Weaver, Carnegie Inst. Technology, Pittsburgh, Pa., 1949. 
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the sill. The force then slowly increased constant value the jump was 
moved further upstream. was suspected that the change force the 
sill was due change the velocity distribution from one end the jump 
the other illustrated Fig. 19. Observations the action glass 
sided flume partly confirmed this theory. 

The nonuniform velocity distribution increased the momentum above that 
existing for uniform velocity distribution the same tailwater depth. 
With the increase momentum the pressure requirement was reduced, thereby 
decreasing the depth. This condition may verified substitution the 
pressure-plus-momentum equation which must satisfied. 


Intial depth 


Force Sill 


Position Jump with Respect Sill 


analysis the results showed that: (a) These tests proved that placing 
sill the end hydraulic jump reduced the tailwater depth below the 
theoretical depth required form jump; and the reduction tailwater 
depth was caused change the velocity distribution due the introduc- 
tion the sillinto The results further disclosed that the assumption 
uniform velocity such was made Messrs. Forster and Skrinde was not 
entirely justified. 

The tests conducted Carnegie Institute were per- 
formed under conditions similar the tests made the State University 
Iowa Iowa City, reported this paper, with one exception. the 
former attempt was made maintain the downstream depth and the tail- 
water the jump the same level instead using the sill raise the tailwater. 
This attempt was not always successful, but was never unsuccessful that 
the sill became the control. The control was located well downstream from 
the sill all times. 

view the fact that nonuniform velocity distribution existed the 
tests conducted the State University Iowa, reduction the jump tail- 
water should occur. The theoretical curve should have been determined 
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the expected tailwater and not the tailwater computed the hydraulic 
jump formula. The tailwater that should occur would less than that com- 
puted. Substitution the expected tailwater into the theoretical formula 
will reduce the value thereby lowering the curve nearer the actual 
results. 


checking through this paper the writers discover once that presents 
interesting mathematical exercise. They are the opinion that every studious, 
serious-minded engineer, before applying the results contained the paper 
actual practice, would feel necessary verify the rather unusual and 
difficult algebraic transformations involved the preparation the paper. 
The usual engineer active practice, however proficient may have been 
mathematics while college, rapidly forgets many details, and soon loses the 
facility the professional mathematician. regretful, therefore, that the 
authors the paper did not include few more brief descriptions the steps 
followed its mathematical treatment. 

attempt remedy this deficiency, the writers venture suggest the 
following explanatory notes, particularly the young inexperienced engineer 
who contemplates contributing further original studies this field. should 
remembered that there are, course, many possible variations details. 

the value from Eq. substituted Eq. the resulting equation 
cubic which means that has three values terms the other 
variables. also rather obvious that factor both members 
the equation. Dividing both sides reduces the equation 
quadratic, which may easily solved any the usual methods. must 
remembered that eliminating the factor one the roots the 


Substituting the quadratic becomes 


the use Eq. the second member this equation may also written 
these three possible forms the best one use the 
particular should determined the existing situation. The 
authors apparently chose the third last the three forms for their subsequent 
derivation Eq. 13a. This natural procedure because Eq. expresses 
the mathematical relation necessary for the formation the jump, and 
the value given Eq. may substituted. This substitution intro- 
duces the height necessary for the assumed controlling sharp-crested weir and 
relates the effect the weir the formation the jump. 

following this procedure, however, formidable and 
dering complications arise. squaring both sides Eq. 11, since and 


Cons. Planning Engr. with TVA, Knoxville, Tenn. 
Eng. Aide, TVA, Knoxville, Tenn. 
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expressed 12a and are both binomials, their squares are 
each trinomials and their product contains nine terms corresponding the 
nine terms contained within the large bracket Eq. 13a. These terms 
contain powers and products the two ratios Eqs. 12a and 
essential sooner later eliminate all redundant forms the variables. 


These substitutions result Eq. 13a printed with slight variation the 
third significant places. reduce Eq. 13a standard form necessary 
use the binomial theorem expand the indicated powers, perform the 
indicated multiplications, and remove the denomipators. The expansion 
produces nearly fifty terms, but combining similar terms the number 
reduced eight. The equation then becomes 


plotting curve (a) Fig. represent 41, necessary compute 
few points the method successive approximations. the 
seventh degree both and which means that convenient integral 
value for either variable substituted Eq. there may seven correspond- 
ing values the other variable which will satisfy the equation. Some these 
values, perhaps most them, may imaginary. Fortunately the curve 
appears rather flat that small number points will suffice de- 
termine it. 

convenient begin with points near the origin where the values 
the coordinates will relatively small. Thus, assuming that the 
independent variable and the dependent variable, let Then Eq. 

Solving the quadratic, 1.396 Only positive values are 
probably practical interest. 

resulting seventh-degree equation indicates that there may seven positive 
roots values which satisfy the equation. successive trial found 
that 3.56 one such value. Continuing similarly, sufficient points may 
found complete the drawing curve (a) Fig. 

the writers’ aim emphasize the importance making clear the 
complex mathematical processes used, that the practicing engineer may 
utilize with confidence the curves shown Figs. and 
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DISCUSSIONS 


CONTINUOUS ARCHES AND BENTS ANALYZED 
COLUMN 


Discussion 


tially study the relation between the elastic deformations each member 
comprising it, and the selection sections that will maintain structural con- 
tinuity. the case plane structures, these elastic deformations occur 
response moment, axial thrust, and shear. Except where the member 
short and deep, the elastic deformations resulting from axial thrust and shear 
are generally small and negligible compared those produced bending 
moment. Following this line reasoning, the column analogy, developed 
Hardy Cross, Hon. ASCE, results simple formula for continuity 
moments, which the deformations from axial thrust and shear are neglected. 
The method useful the solution three-redundant structure because 
susceptible systematic computation when arranged tabular form. 

The extension the column analogy the analysis continuous struc- 
tures, presented Mr. Yu, contribution engineering literature 
despite the fact that its range application limited certain cases, the 
same many other methods and modifications the field structural 
analysis. 

its application the analysis arch spans (particularly when the arch 
flat), the column analogy good only for preliminary investigations because 
neglects the effect axial deformation, omission that seems unwarranted. 
other words, unless axial deformation can taken into consideration, 
questionable the author’s conclusion (that his method especially useful 
device for reviewing designed structures) justified the case arch spans. 

Example the author offers continuous bent for analysis. The same 
problem was solved the writer using the method moment distribution. 
With the aid ordinary slide rule, required only hours finish the en- 
tire computation, including the values influence ordinates for moment, 


paper was published February, 1949, Proceedings. Discussion this 
has appeared Proceedings, follows: June, 1949, Jack Benjamin; and October, 1949, 
eslie Meszaros, Eremin, and Polivka. 


Design Engr., Naval Port Facilities, Shanghai, China. 
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thrust, and shear, corresponding the data Although comparison 
with the author’s working time great consequence (the speed com- 
putation depending large measure the degree acquaintance with each 
particular method), the relative accuracy results proper criterion. 
Since presentation other methods beyond the scope discussion, the 
writer’s computations cannot presented; but the results show, invariably, 
smooth curve influence ordinates for each stress function whereas the 
author’s results show marked divergence the remotest spans. The errors 
are most likely mistakes computation because systematic errors com- 
putation would produce influence ordinates that revealed general disagree- 
ment, rather than the isolated ones actually produced. clear that the 
column analogy limited basis for analysis because offers method 
applying independent check. Mistakes computation will carried 
along and thrown into the final result unless (as demonstrated the author) 
they can discovered the equality reciprocal deflections, which 
check not available until the computation completed. For the case 
long series continuous spans, the successive combination spans necessary 
arrive single span wherein the particular stress functions are studied 
would entail many tabulations that slight mistake computation, not 
rigidly checked the end each process combination, may result 
final, cumulative, gross error. 

Furthermore, getting influence lines using the Mueller-Breslau prin- 
ciple, the computation deflection ordinates remote spans, due unit 
strains corresponding the desired influence functions, permits rather high 
degree precision shown the significant data Fig. the writer’s 
solution the same problem, uniform degree precision was controlled 
using ordinary slide rule. tracing influence lines for stress functions 
sections other than those illustrated the author, necessary find 
many independent sets indeterminate moments column analogy (moments 
due applied unit strains particular section) there are numbers 
sections intersected, and then compute the deflections. moment dis- 
tribution, only matter simple arithmetic applied the principle 
statics, compute influence ordinates any number sections, using only 
the influence lines for fixed-end moments and the joint moments distributed 
from joint with arbitrary unbalanced moment. 

The column analogy computation gives the combined effects sidesway 
and joint rotation but affords clue the proportions between them. 
Moment distribution yields that result directly and much preferred. 

Example the paper shows clearly the merit Mr. Yu’s method. The 
directness with which the substitute members are found simplifies the solution. 
The same problem, solved the method moment distribution, would 
require two sidesway corrections and two simultaneous equations, making 
perhaps less favorable than the author’s method. 

summarize, the writer feels that, although the author’s contribution 
engineering literature highly commendable, its field application the 
analysis continuous structures limited, being neither rapid nor 
accurate many other methods. 
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DISCUSSIONS 


WIND STRESSES SEMI-RIGID CONNECTIONS 
STEEL FRAMEWORK 


Discussion 


GEORGE MURPHY 


Assoc. ASCE.—The determination the actual 
strength semi-rigid connections depends (as the author states) the avail- 
ability experimentally determined empirical curves, one for each joint 
combination, assuming identical quality workmanship amd materials 
every case. Therefore, the problem discover what limiting assumptions 
will keep one the safe side uncertainty. 

Seldom does ordinary design procedure beyond computations for the ten- 
sion, bearing, and shear the rivets welds, while assuming the base material 
the members and the connection pieces 
fully nondeformable—despite full awareness, 
the designer, the flexibility and plasticity 
the paris standard connection. The frame 
(excluding wind bents) designed assuming zero 
end moments (that is, assuming that there 
restraint against the pulling away the top 
connection); and yet, subsequently, the super- 
imposed concrete floor designed the framing 
connections were fully continuous. Running 
cracks the concrete floor, along the steel girders 
which the steel floor beams are framed, are not 
unknown occurrences. They bear witness the 
ill effects this contrasting approach the de- ror 
sign different parts the same structure; CuRVE 
any effort correct this approach should en- TIME 

EVERSAL MOMENT 
couraged. 

The general type stress diagram for connections (Figs. and 6), which 
moments are plotted against angular deflections, might questioned the 
area positive negative increase from zero moment, demonstrated 


paper Basil Sourochnikoff was published February, 1949, Proceedings. 
this paper has appeared Proceedings, follows: September, 1949, Batho. 
With Samborn, Steketee, Associates, Toledo, Ohio. 
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Fig. 13. analogous illustration would the free play the steering 
column car boat when the direction the turning force reversed. 
break the elastic curve, such this, would have considerable effect 
altering the distribution moments the structure where not all the joints 
would subjected reversal moment. Stress states and would 
then not identical. 

defining the conditions for keeping connection within the elastic range, 
the authof’s association algebraic and absolute numbers may lead con- 
fusion. This danger particularly apparent when one considers that given 
connection may have different elastic limits Fig. 13) for different direc- 
tions rotation. Eq. 12a would more mathematically precise the form: 


(22) 


other words, algebraically greater than the negative the absolute 
value the elastic limit the joint when rotated the negative direction. 

Conclusions (1) and (2) are undeniably correct. Conclusion (3), however, 
qualitative conclusion within more less indeterminate limits; and 
caution required avoid its improper application specific cases. Judg- 
ment required estimating the capacity inherent the connection 
undergoing inelastic deformations without failure, and estimating 
the magnitude the inelastic deformations will required undergo. 
Furthermore, the statement that “adequacy” the connection for combined 
gravity and wind loads will found resist wind moments 
alone found, when the latter are computed assuming elasticity the con- 
nections, contingent, not only upon sufficient deformability the connection, 
but also upon sufficiently large ratio the wind moments those caused 
the gravity loads. Otherwise, windless climes, connections would 
required. This recognized the design specification that, any case, the 
connection for the gravity loads. 
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EFFECT DRIVING PILES INTO SOFT CLAY 


Discussion 


HOUSEL 


this paper, was conducted the Union Station Detroit, Mich., connection 
with the design grade separation structures the John Lodge Express- 
pertinent compare the conclusions drawn that case with 
those presented the authors. some respects the later investigation con- 
firms the authors’ deductions and provides some additional data regarding 
this problem about which experimental results are still rather fragmentary. 


The primary objective the tests was study methods controlling the 
volume displacement due pile driving which sometimes has harmful effects 
adjacent structures utilities. With respect this phase the investiga- 
tion was found that the use tapered, open end steel piles the type used 
did eliminate major part the volume displacement. Subsequent loading 
tests indicated that those piles which the clay core was not removed had 
ultimate carrying capacity controlled the soil comparable that mea- 
sured the piles filled with concrete. 

Largely result the publication the earlier investigation made 
the authors was decided conduct some parallel studies connection with 
the Union Station pile tests study remolding and moisture content changes 
due pile driving. The tentative conclusions drawn from the investigation, 
and which may compared with the authors’ deductions, are follows: 


The Detroit clay deposit consists least two depositions lacustrine 
clay showing different reaction pile driving far loss shearing 
resistance concerned. The physical characteristics these deposits, de- 
termined index tests, show some similarity those reported the authors. 


1948, Proceedings. Discussion this paper has appeared Proceedings, follows: May, 1949, 
Leonardo Zeevaert, Jacob Feld, Joseph Fegan, and Gregory Tschebotarioff; June, 1949, Philip 
Rutledge, Skempton, James Schuyler, and Kantey; September, 1949, Krynine, 
and Robert Leggett; October, 1949, Stuart Avery, Jr., and Stanley Wilson, and George 
Freeman and James Parsons; and November, 1949, LeRoy Crandall. 

Associate Prof., Civ. Eng., Univ. Michigan, Ann Arbor, Mich. 


Burkey, Proceedings, International Conference Soil Mechanics and Foundation Eng., Rotter- 
dam, 1948, Vol. pp. 146-154. 
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Above approximate El. 530 the Union Station, clay within in. the 
pile suffered loss shearing resistance comparable the value when remolded. 
Determination recovery was based two types shear tests made approxi- 
mately 2.5 months after pile driving. Below El. 530 there was apparent 
loss strength due pile driving distance approximately in. from the 
pile. tests were made determine recovery the soil with time. Shear 
values the lower clay deposits 2.5 months after driving were substantially 
higher than the original shear strength found the check boring but within 
the range variation found series exploratory borings. These results 
were confirmed two similar borings distance in. from the piles. 

two borings distance several feet from the piles there was 
evidence loss shearing resistance due remolding either the upper 
the lower clay deposit. one case shear values the lower deposit were 
greater than the check boring but within the range found the series ex- 
ploratory borings. This would tend confirm that the higher shear values 
found the lower clay were the result natural variation the deposit 
rather than alteration produced pile driving. 

There was positive evidence any change moisture content 
either the top bottom clay deposit but one depth where there was 
particularly soft and highly plastic clay there were scattered results from which 
might concluded that decrease moisture content had taken place 
within in. the piles. 


presenting these data the writer wishes point out that, although the 
results reported are based duplicate borings and duplicate laboratory shear 
tests and while they indicate some trends, there such wide variation 
shearing resistance and physical characteristics the natural state that the 
statistical basis comparison becomes extremely narrow, making the available 
data comparatively fragmentary. For example, both investigations there 
are only data two points with respect distance from the piles, one ap- 
proximately in. and one several feet, giving little measure the boundary 
zones within which remolding important. Perhaps more important the 
possible variation with time, particularly the rate recovery the loss 
strength. this point there quite incomplete coverage the results re- 
ported from both investigations, and considering the character the tests and 
the expense involved not likely that such information will made avail- 
able under ordinary circumstances. 

would appear the writer that there are other, and perhaps more im- 
portant, aspects the effect driving piles soft clay which attention 
might directed and which observation and experience may provide more 
conclusive information. For one thing, matter common experience 
pile driving that the load-carrying capacity the piles determined 
conventional pile driving formulas may very small the time the driving 
completed, but within period hours few days friction piles 
and offer much greater resistance penetration. fact, common 
practice take advantage this well-known phenomenon meeting specifica- 
tions that require certain pile capacity determined commonly used 
formulas. 
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Some years ago the Detroit Edison Company conducted rather extensive 
series load tests friction piles determine the recovery the soft clay 
from the detrimental effects pile driving. These tests were made direct 
loading two 16-in. (outside diameter) steel piles, approximately long, 
which were driven the company’s Delray Power Plant Detroit. Soil con- 
ditions the area are comparable those under discussion the authors, 
being the River Rouge area Detroit where the entire depth clay 
hardpan consists soft-plastic lacustrine deposit. 

Immediately after the piles were driven, load testing equipment was set 
and direct loading tests were conducted intervals starting from about 
hours after the piles were driven and continuing for about days. The 
results these tests have never been published but seem pertinent the 
subject now under discussion and provide evidence the nature the 
temporary loss strength and subsequent recovery such clays after pile 
driving. Observations these phenomena indicate rapid recovery 
resistance accompanied parallel changes elastic properties 
which are seldom associated with changes moisture content and consolidation 
the soil mass. 


Load in, Thousands Pounds 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Inches 


Test No. 20, Micu. 


typical load settlement diagram taken from the series thirty-six tests 
conducted the two piles shown Fig. 35. This diagram indicates the 
resistance characteristics the clay mass measured loading so-called 
friction pile. Even though the clay soft highly plastic gives evidence 
extended range elastic behavior shown the linear relation between 
load and settlement extending the load limit. After the load limit has 
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been exceeded there stage rapid progressive settlement shown 
almost horizontal load settlement line followed rebound elastic recovery 
after the release the load. Throughout the series tests was found that 
the rebound was approximately equal the total elastic deformation deter- 
mined the intersection the elastic load settlement line with the failure 
line, point which designated the load limit, Pz. The shearing limit, Ps, 
has been selected the load carried deformation equal the total elastic 
deformation. The difference between the shearing limit and the load limit 
the curved portion the load settlement diagram considered the point 
resistance. this connection should pointed out that the friction pile 
Fig. was tested with the clay core inside the pipe extending almost the 
top the pile, condition resulting immediately after driving. Fig. 
shown comparison between rebound and total elastic deformation for the 
series tests such friction pile, supporting the conclusion that the re- 
bound and total elastic deformation were equal. 

Fig. shown load settlement diagram for the same pile after the 
clay core had been removed for the purpose completely eliminating point 
resistance. the first place evident that the difference between the 
shearing limit and the load limit has almost disappeared and the failure the 
pile quite abrupt. greater importance, however, the comparison 
afforded immediate repetitions load shown Fig. type test 
was made number times throughout the series tests demonstrating 
beyond any reasonable doubt that friction piles such cohesive clays could 
loaded rapid progressive settlement without impairing the shearing resis- 
tance the material. 

The variation load-carrying capacity the pile, which may assumed 
equivalent the shearing resistance the soil, provides the most im- 
portant data available from these tests. This variation shown Fig. 
which the shear limit previously defined has been plotted for one the 
piles with the time days since driving the abscissa. few hours after 
driving the shear limit was approximately 25,000 lb. This increased rapidly 
and period days the shear limit had increased approximately 
77,000 lb, about three times its initial value. this point the pile had 
developed capacity approximately equal that which had been estimated 
from the independent shearing resistance tests which the average shear for 
this depth was approximately 275 per ft. duplicate test the second 
pile this time provided very satisfactory confirmation the total shearing 
capacity. 

One the most interesting observations which this series tests provided 
resulted from what might termed accident that occurred after the first 
days had elapsed. the preceding piles was tested rapid pro- 
gressive failure and this point became necessary put additional 
length pipe this pile order continue the tests. The construction 
superintendent who had been assisting the project was very much interested 
and very anxious everything possible bring the investigation 
successful conclusion. Consequently, order eliminate any possible settle- 
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ment from incomplete seating the new length pipe the joint, brought 
the pile driver back the job and drove the pipe about in., possibly striking 
more than six blows. The surprising result shown the curve Fig. 38. 
The next test after redriving the pile shows abrupt decrease shearing 
resistance which was substantial although not great the original loss 
strength. 
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Time, Days 


seems obvious that the slight movement the pile involved redriving 
could not have produced any significant amount remolding the soil 
change the moisture content. should also pointed out that the dis- 
placement involved was only little more than that which took place each 
the static load tests which the soil was carried rapid progressive dis- 
placement. may stated with considerable certainty that deforming the 
soil (by rapid shearing displacement under static load) did not impair its 
strength; but the shock just few blows the hammer produced marked 
decrease strength. Furthermore, the decrease shearing resistance was 
accompanied parallel decrease elastic deformation recorded both 
the rebound and the total elastic deformation. These data are shown Fig. 
and, the writer’s opinion, present very strong evidence that the primary 
source the loss strength and its subsequent recovery can traced back 
the very source cohesion adsorption molecular orientation soil mois- 
ture. The phenomena may called thixotropic effect, which character- 
istic colloidal systems and has been observed clays the character under 
discussion. 

This series tests produced evidence least one other interesting 
phenomenon. may mentioned passing even though its basis some- 
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what speculative. the tests were conducted the shearing resistance in- 
creased with time maximum days shown test No. (Fig. 
38) where the shearing limit approximately 90,000 lb. This total capacity 
was substantially greater than had been anticipated from previous shearing 
resistance tests and greater than the shearing resistance the duplicate pile 
that was being tested. the same time the test data gave evidence 


Pipe driven 


Core removed from pipe 


0.20 


Time, Days 


Rebound and Elastic Inches 


disparity between elastic deformation and shearing resistance which led 
the speculation that dynamic point resistance was being built result 
repeated tests which produced rapid progressive penetration the pile. The 
magnitude this dynamic point resistance was quite high, approaching the 
passive pressures that might anticipated under the static head some ft. 

was then decided two things obtain further evidence this 
phenomenon. the first place, the pile was left without load repetition for 
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approximately days and then tested shown test No. (Fig. 
total shearing capacity that time had dropped back approximately 77,000 
This was fairly close that which had been anticipated from other test 
data, indicating that the high intensity pressure developed the point 
the pile had been dissipated. the succeeding test No. 28, however, this 
dynamic resistance was substantially recovered. the next test the core 
was removed from the pipe order eliminate point resistance completely. 
This removal was accompanied considerable disturbance caused the 
inflow soil into the bottom the pipe, was found the soft clay that 
either head water substantial core was required prevent such loss 
ground. the subsequent tests, shown Fig. 38, the shearing capacity 
the pile was recovered and the series tests was terminated. 

Acknowledgment.—The tests the Union Station, referred this dis- 
cussion, were conducted jointly the Michigan State Highway Department, 
the Wayne County Road Commission, and the Union Metal Manufacturing 
Company. 
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MULTIPLE-PURPOSE RESERVOIRS 
SYMPOSIUM 


Discussion 
DANA AND CLARENCE JARVIS 


Dana ASCE.—Considering the papers this Symposium 
whole, the reader immediatley impressed the number and authority 
agencies having interest federal multiple-purpose development 
particular, and nearly equal interest private developments. Assuming, 
for example, that the Corps Engineers, United States Department the 
Army, under Congressional authorization, has paramount interest planning 
and design, because flood control navigation functions both, must: 
(1) Obtain the approval the Federal Power Commission the design and 
the comprehensiveness the river development scheme; (2) cooperate with 
the Bureau irrigation involved arrange for the United 
States Department the Interior sell any power available; (3) perhaps 
coordinate power planning with private utilities independent market 
not available; (4) secure recommendations from the Fish and Wildlife 
Service the Department the Interior; (5) determine the extent participa- 
tion the National Park Service the same department with respect rec- 
reational possibilities; (6) procure the hydrometeorological services the 
United States Weather Bureau and the United States Geological Survey; (7) 
ascertain public health requirements; and then (8) receive cooperation and 
assistance lower state and municipal levels. Historic and archeological 
investigations must included the National Park Service, assisted the 
Smithsonian Institution Washington, C., and perhaps local uni- 
versities. The soil conservationists have interest. this requires special 
appropriations for each agency, shifting funds from one the other, 
necessitating government budget requests and reviews, special agreements. 
How complicated the life the government planning and design engineer 
becoming and how divided the authority! 


Symposium was published March, 1949, Proceedings. Discussion this paper has 
appeared Proceedings, follows: June, 1949, Herman Hilton, Wesley Gilbertson, Roger 
Amidon, Bishop, and Davison; and September, 1949, John Cotton, and Abraham 

Chf., Power Studies Branch, Div. Water Control Planning, TVA, Knoxville, Tenn. 
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Little has been said about the necessity obtaining local approval, 
coordinating, and compromising the appeals special pressure groups, 
overcoming selfish and individual attitudes all levels. course, the 
greatest good the largest number individuals supposed the con- 
trolling theory, but political and financial controls can violate this. One the 
greatest evils today unwillingness compromise issues any level 
from the individual even the international. The question too apt 
raised, ‘‘How does this affect the individual, the community, the 

Yes, national water policy needed, but how one obtained? 
Congress ready modify its stipulations that one agency sell its power 
the switchboard while another builds its transmission lines the market, 
any particular agency willing depart from its own pet theory alloca- 
tion costs multiple-purpose objectives determining project economics? 
flood control greater less importance than power except for reasons 
political expediency? Are local benefits sufficiently large outweigh local 
taxes government projects? Should theoretical rather than actual interest 
charged against government project which competes with private business? 
Will nonreimbursable appropriations made the expense the United 
States large that purchasers water and power will not required 
pay for services provided for the general public the same restricted area? 
Shall large landowners receive the same benefits from irrigation small 
landowners? Shall residents within flood control district contribute the 
project some equitable basis shall the federal government bear the entire 
expense? Shall governmental planning allocated numerous bureaus 
consolidated new bureau? One can indefinitely raising ques- 
tions requiring compromise which probably never will actually compromised 
until the character the individual and the groups improved and selfish 
interest largely eliminated. 

national water policy must flexible because two cases are exactly 
alike. one case, may necessary sell power the switchboard 
the one available utility, thereby continuing monopoly, because any other 
market could not built readily. another case, the public interest 
might hurt following such policy. The residents within proposed 
flood control district might not have the financial ability participate, 
matter how desirable the flood protection might be, because already existing 
financial responsibilities. Many municipalities are overburdened with heavy 
debt, perhaps result mistaken past policies, and cannot undertake badly 
needed improvements. This condition, course, raises the important question 
whether community which has been thrifty maintaining 
you-go basis should now penalized assist through taxes less thoughtful 
community. Water laws and practices one section the United States 
differ from those another. Basic fundamental questions are woven all 
through the problem national water policy. 

Desirable some national water policy may be, seems likely that will 
long time before one can developed. the meantime, the step-by-step 
progress made recent years should continued. Indeed, any fixed policy 
may undesirable because the constantly changing conditions the 
United States. 
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ASCE.—Relatively few water storage reservoir 
projects large capacity and major importance have been developed serve 
single purpose. Ordinarily one principal interest was stressed, with inci- 
dental benefits others, and possibly some disadvantages and inconvenience 
accrued unless prior rights were recognized and fully protected. 

The legislation embodied House Document No. 308,** which was enacted 
into law Section the River and Harbor Act January 21, 1927, initiated 
new era practical hydrology. directed the Chief Engineers prepare 
technical reports the principal river systems the United States, except 
the Colorado River (on which the Bureau Reclamation, United States 
Department the Interior, had already made notable progress), from the 
fourfold standpoints navigation, flood control, power development, and 
irrigation. Naturally, the interests municipal water supply, manufacturing 
industries, recreation, fish and game life, and other outstanding items 
national life and economic structure were duly investigated and included the 
final reports. The impressive array Document No. 308 Reports,”’ 
published during the following decade, helped establish the concept 
“multiple-purpose which both the principal and the incidental 
beneficiaries were recognized, along with their eligibility and probable ability 
share the capital cost the undertakings. This new concept has made 
possible the launching several important key projects the realms 
control, navigation, hydroelectric power development, irrigation, and municipal 
water supply. Without the hydroelectric energy made available such 
projects Grand Coulee Dam (on the Columbia River), Bonneville Dam (on 
the Columbia River), Hoover Dam (on the Colorado River), the notable 
features the Tennessee Valley Authority, and other monumental achieve- 
ments, just how could this nation have accomplished its tasks during the 
years World War II? 

substantial part the total costs such projects (with their demon- 
strated strategic value and importance during national crises) must charged 
the general welfare, while subscribers under irrigation, municipal water 
supply, and hydroelectric power projects pay all that can rightfully assessed 
against them, such basis financing both justifiable and commendable. 

Alleged conflicts major interests multiple-purpose reservoir project 
are too often grossly overdrawn. Conceding that the flood control ideal 
inclines toward empty reservoir certain seasons, and that hydroelectric 
power interests prefer high storage head, both ideals must readily submit 
limitations and modifications their own behalf, and should make even 
greater concessions the spirit compromise promote conservation and 
utilization basic natural resources. properly designed and operated, 
reservoir filled nearly spillway level proof mission accomplished for 
that particular season period. still capable storing and automatically 
releasing regulated volume floodwater over the spillway, with but minor 
negligible damage the pasture woodland below the designed high-water 
contour, for which flowage rights are already secured. remarkably large 


Cons. Engr., Salt Lake City, Utah. 
House Document No. 308 69th Cong., ist Session, 1928. 
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proportion the total reservoir capacity ordinarily located above the spill- 
way level. With assurance that automatic release will continue the prede- 
termined rates, over the open spillway, the operators are relieved needless 
worries and the pitfalls immature decisions regarding emergency drawdown. 

Recent advances the field hydrometeorological services the United 
States Weather Bureau and cooperating agencies have effectively removed 
the probability, not the possibility, unwarned approach catastrophic 
storms (and resulting floods) except those associated with structural failures. 
This situation certainly justifies the adoption intelligent and liberal policy 
the field reservoir drawdown. With long-range weather forecasts attaining 
higher percentages accuracy and dependability, becoming more justi- 
fiable maintain considerable storage volume floodwater for orderly release 
promote beneficial use. the other hand, demand for hydroelectric energy 
may justify drawdown below the ideal storage level when there reasonable 
assurance increased volume river discharge the near future, capable 
replenishing the storage the desired level even above that level. 

Some the most important and valuable storage capacity reservoir 
may consist what known bank storage. Lateral infiltration into coarse 
granular deposits other readily permeable geological formations above 
Elephant Butte Dam the Rio Grande New Mexico accounted for return 
flow, during the final drought year series ending 1940, which was suffi- 
cient offset all reservoir losses and yield additional dividend approxi- 
mately 47,800 acre-ft. The 46-year mean annual loss had been calculated 
198,700 acre-ft from records the Bureau Reclamation, summarized, 
tabular form* and graphical form,* the writer 1943. some cases, 
geological reconnaissance may determine the advisability opening trenches, 
cuts, shafts penetrate the potential underground storage strata, especially 
where the topsoil mantle would retard the inflow and outflow unduly. 

The Symposium comprehensive scope and constitutes valuable con- 
tribution the subject multiple-purpose reservoirs. 


Gleanings from the Field Hydrology,” Clarence Jarvis, Brentano’s Inc., 
1943, pp. 40-41. 
37. 
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DISCUSSIONS 


ENGINEERING CONCEPT FLOW PIPES 


Discussion 


ANDRE JORISSEN, WESTON GAVETT, AND 
ARTHUR COLLINS 


worthwhile endeavor approach the problem 
flow pipes natural roughness synthetically made Mr. Harris. 
The main obstacle comprehensive solution lies the complexity the 
roughness pattern and the resulting difficulty defining, and hence 
measuring, roughness. 

attempt made characterize roughness specific roughness 
coefficient. should understood that this coefficient the constant 
friction factor fully turbulent flow pipe unit diameter, for given 
roughness pattern. specification the unit length used necessary 
stated that the value Eq. must expressed the same units. 

The analysis transition from laminar fully turbulent flow based 
derived from the author’s experimental study smooth brass pipe. 
This formula similar relation which was propounded, early 
1922, Charles Liege, Belgium: 

0.532 

The Hanocq formula has been widely and satisfactorily used western 
Europe since the early 1920’s. based numerous experiments, both 
with air and water and carried the Reynolds numbers 500,000. Also, 
has the advantage pure cube root the computation the second term. 

Reynolds number 3,240,000 having been reached the Nikuradse 
experiments, desirable investigate what extent classical formulas 
agree with experimental results. Typical formulas for smooth pipes (in 
addition the Hanocq formula already mentioned) are follows: 


paper Charles Harris was published May, 1949, Proceedings. Discussion 
this paper has appeared Proceedings, follows: September, 1949, Harrington; October, 1949, 
Ralph Powell, Stevens, George Whetstone, Hickox, and James Halsey and William 
Owen; and November, 1949, Tsen-Ding Chen, and Lynn Perry. 

Associate Prof. Civ. Eng., Pennsylvania State College, State College, Pa. 


propos calcul des pertes charge dans les conduites Charles Revue 
Universelle des Mines, Liége, Belgium, February, 1922, Series pp. 217-234. 
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1.14+2 log 


Values 


LEGEND 


Nikuradse 
Jorissen 


1.2 1.4 1.6 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 


Values Log 


Fia. 


The formula (1913)— 


for 100,000; the Lees (1915), based experiments 


0.6104 


for 460,000; and the formula (1930)— 


(27) 


dem Gebiete der No. 131, 

the Flow Viscous Fluids Through Smooth Circular Charles Lees, Proceedings, 
Royal Soe. London, Series Vol. 91, 1915, 46. 

Pannell, Philosophical Transactions, Royal Soc. London, Series Vol. 214, 1914, 199. 

von Platte und bei hohen Reynoldsschen Schiller and Hermann, 
Ingenieur Vol. 1930, 392. 
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—which gives results good agreement with the Nikuradse experiments 
Reynolds number 2,500,000. discrepancy however apparent for 
50,000. The formula (1932)— 


valid for Reynolds number greater than 100,000, which the limit 
validity the Blasius formula. 


TABLE FORMULAS FOR FRICTION Loss 


Reynolds’ NIKURADSE Harris 
No. (R) HERMANN 


(writer) 


Eq. Eq. Eq. Eq. Eq. 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 


More recently, the writer has derived comprehensive formula based the 
Nikuradse measurements velocity distribution and pressure loss. This 
expressed 


and checks remarkably well with the experiments Nikuradse and the writer 
for all values the Reynolds number. 

comparison the various formulas (see Table shows fairly good 
agreement the range validity each formula, but the agreement not 
extended over this range. and give close values the high 


der turbulenten Strémung glatten Nikuradse, Forschungsheft, 
Verein Deutscher Ingenieure, No. 356, 1932. 
des pertes charge continues dans les conduites Andre 


Jorissen, Bulletin Centre d’Etudes, Recherches d’Essais Scientifiques des Constructions Génie 
Civil Fluviale, Liége, Belgium, 1949. 
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2.5 
LEGEND 
12% 
pipe numbers 


log 


Values 


0.5 


R. 
Values log 


Reynolds numbers, slight discrepancy being noticeable the extreme values 
the latter parameter. However, for all practical purposes, the two formulas 
may said agree. 

The proposed method characterizing roughness based the con- 
sideration the specific roughness coefficient. doubtful, however, such 
coefficient will prove satisfactory the study commercial rough pipes. 
the experience the writer, the mean height roughness protuberances 
may vary considerably from one pipe another, even for new pipes ap- 
proximately equal diameters manufactured the same factory. This opinion 
was obtained measuring the mean roughness height pneumatic 
was shown that, the new industrial steel pipes studied, 
characteristic value for given pipe. The value may vary from one 
section another and from point point one section but, averaging 
large number measurements, characteristic mean value obtained. 
stated before, differs considerably from one pipe another. This stresses 
the difficulty using any the classical formulas for rough pipes and seems 
restrict the use the specific roughness coefficient. Pneumatic feeler 
measurements, the contrary, are easier make than pressure loss determi- 
nations and may further repeated from time time follow the evolution 
the pipe roughness. 

Another conclusion the aforementioned study was that, the zone 
fully developed turbulence, the value was given with very satisfactory 
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values log 


and White 
(homogeneous roughness) 


e 
Values 


approximation the Nikuradse formula: 
(30) 


which was the mean roughness height, measured the pneumatic 
feeler (Fig. 22). 

smooth transition the law variation function the Reynolds 
number was recorded. This also agreement with the Colebrook and 
White experiments. this zone, the law variation for industrial pipes 
situated between the Nikuradse law for homogeneous roughness and the Cole- 
brook and White law for heterogeneous roughness. The writer’s experiments 
show classification following the roughness homogeneity (Fig. 23). 


Gardner Williams and Allen Hazen, Members, ASCE, has been described 
the author (under the heading, Some Early 
“based fallacy that should have been apparent the time.” Mr. Harris 
further remarks that: 


formula was developed, selecting only the data 
from the smoothest pipes attainable, the exponent reducing 1.85. 


was then assumed that all rough pipes would follow that same exponential 
aw’ *x* * 


Associate, Clyde Potts, New York, 


Gardner Williams and Allen Hazen, John Wiley Sons, Inc., New 
York, Y., Ed., 1920. 


pipe Nikuradse (homogeneous roughness) 
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the Messrs. Williams and Hazen 
pointed out that— 


the values the exponents vary with different surfaces and 
also their values may not exactly the same for large diameters and for 
small ones, nor for steep slopes and for flat ones. 

“The exponents the formula used were selected representing 
nearly possible average conditions, deduced from the best available 
records experiments upon the flow water such pipes and channels 
most frequently occur water-works practice.” 


That the Hazen-Williams formula fact represents fair average for pipes 
water systems indicated its continued use and the results many 
tests before and after its publication. The convenience such exponential 
formula for solving network problems, using the ‘‘Hydraulic Tables” and the 
Hazen- Williams slide rule generally recognized. 

The exponent 1.85, between the usual values 1.75 for smooth pipe 
and 2.00 for complete turbulence, seems represent most probable value. 
With the increasing use smoother pipe surfaces, the exponent 1.85 may 
become high rather than low. 

For important transmission mains, particularly those large diameter, 
the water works engineer will study all available data and theory; but, for 
distribution system studies, the Hazen-Williams formula should continue 
favorite tool. 


Assoc. ASCE.—The writer’s concept the 
flow water pipes best illustrated multitude points logarithmic 
chart. They are numerous that they outline band. Each point has 
been established two readings, velocity and loss head, reported 
observer. 

Each new formula proposed for pipe flow the result laboratory test 
definitely fixed line within the chart band. There can about 
many different formula exponents there are stars the Milky Way. Re- 
gardless the formula used the data are someplace the chart and choice 
can made after reading the testimonials. 

Other factors enter into flow measurements and characteristics the field 
which are not present the laboratory where the setup includes small pipe 
with perfect bore. Undulations and ovalness (out round) are conditions 
that can have marked effect upon energy losses. 

The study the velocity distribution velocity factor should in- 
cluded pipe experiments method classifying pipe conditions and 
making the data more comparable. reasonable assume that some 
time the pipe line will designed provide certain velocity factor which 
will indirectly insure maximum transportation efficiency. 


Hydraulic Tables,” Gardner Williams and Allen Hazen, John Wiley Sons, Inc., New 
York, Y., Ed., 1920, 


Cons. Engr., Berkeley, Calif. 
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DISCUSSIONS 


MINERAL LAND SURVEYS 


Discussion 
ARTHUR MCNAIR 


Assoc. ASCE.—In this paper the writer pre- 
sented material the history and methods making mineral land surveys 
because felt that the subject had not received any attention from the engi- 
neering profession since the turn the twentieth century. The comments 
emphasize the timeliness the subject because the old-timers familiar with the 
procedures peculiar this type surveying are rapidly diminishing numbers. 
Furthermore, the number both patented and unpatented mining claims has 
been increasing continually during these years while the same time the 
corner monuments fixing their positions the ground have become increasingly 
difficult find. 

This situation, stated Mr. Meldrum, creates fertile field for the 
development young engineers. large sections the old mining districts 
are consolidated for unified operation larger companies, the need for resurvey 
the old claims, together with the location and patenting new claims 
cover the vacant fractions, immediately felt. more than passing 
interest note that titles claims old established mining district are 
usually such condition that least two three years work required 
legal counsel before mining company can consolidate ownerships and begin 
operations. 

The discussions indicate that there are continual changes the administra- 
tion the lands the public domain and that these have produced some 
changes the procedure making mineral land surveys. substantial 
contribution this paper has been made the discussions which describe the 
functions the Bureau Land Management, United States Department 
the Interior, which now administers the lands the public domain. 

Because the discovery gold—and other mineral values—was frequently 
the reason for the early settlement many the western states, and because 
each mining district grew autonomously, not surprising that today 


paper Arthur McNair was published November, 1948, Proceedings. Discussion 
this paper has appeared Proceedings, follows: April, 1949, Charles Eager, and Robert Lenon; 
and June, 1949, John Meldrum, and Arthur Kidder. 


Associate Prof. Civ. Eng. and Head Surveying Dept., Cornell Univ., Ithaca, 


1539 


ch 


1540 McNAIR MINERAL SURVEYS Discussions 


there are local variations mining laws and regulations. The hope the 
writer that these variations would presented those discussing the paper, 
thereby making more applicable the entire United States, has been 
least partly fulfilled. Since California probably does more placer mining 
than any other state, the discussions presenting the variations placer surveys 
that state compared others are important. interesting 
that the possessory right unpatented claims the form locations tax- 
able California this sets different economic risks for locator than are 
established where locations are not taxed. 

The requirement Arizona that end-center monuments (or their 
counterparts each end the lode line when the lode does not cut the end 
line exactly its center) must set differs from the procedure 
other states. The writer cannot understand the necessity for end centers. 
Certainly they serve useful purpose, and one cannot help but surmise 
what the mineral surveyor might designate the lode line for claim located 
region disseminated copper porphyries where there are precise 
boundaries, such Arizona, region sedimentary so-called blanket 
deposits such occur many places the western states. 

With the introduction European-make transits into the United States, 
considerable part the writer’s discussion concerning the precision 
mineral land surveys will modified. Some these instruments are not 
only much lighter and easier carry but are more accurate and much easier 


read than the conventional American transits described the paper. 
phase surveying are portability, accuracy, and ease reading more essential 
than mineral surveying, and the next few years will doubtless see im- 
provement procedures corresponding the improvement instruments. 
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DISCUSSIONS 


DIFFUSION SUBMERGED JETS 


Discussion 


HAROLD 


Jun. ASCE.—The value this paper lies not only 
the specific results obtained, but the method attack, which may 
applied with modification necessary submerged jets different nature 
and application. For instance, one the specific examples referred the 
authors, efflux from submerged sluice gate, cited for comparison. The 
upper limit submergence sluice gate (that is, infinite vertical extent 
tailwater) may represented substituting boundary for the plane 
symmetry the completely submerged two-dimensional jet shown Figs. 


Energy; 


turbulence 


d 


Fic. REPRESENTATION FLOW FROM SuBMERGED 


paper Albertson, Dai, Jensen, and Hunter Rouse was published 
December, 1948, Proceedings. Discussion this paper has appeared Proceedings, follows: June, 
1949, James Holdhusen, Duilio Citrini, and Stanley Corrsin; and September, 1949, Douglas 
Baines, and Abraham Streiff. 


Instr., Civ. Eng. Dept., Georgia Inst. Technology, Atlanta, Ga. 
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and 21, whereas the lower limit represented free hydraulic jump 
forming this boundary downstream from the efflux section. far the 
majority cases involving flow under sluice gate occur nearer the limiting 
case the free jump: Fig. represents, schematically, typical intermediate 
condition. 

The authors’ investigation serves indicate the flow pattern for conditions 
approaching the upper limit. the other hand, lower degrees sub- 
mergence, the proximity the free surface cannot ignored. The new 
effects introduced are: (a) definite gravitational effect evidenced 
the depression the surface indicated Fig. 31, and which may considered 
the analysis the inclusion Froude number; and the limiting the 
fluid downstream from the efflux section finite vertical extent which calls 


H: 


Distribution Intensity Turbulence 


for variable describing the relative depth the tailwater. The foregoing 
opening the tailwater depth; the mean flow per unit width 
gate opening; and the acceleration gravity. Hence, when modified 
apply two-dimensional flow from submerged sluice gate, the formula for 
relative velocity any point the flow (Eq. would take the form: 


which the mean velocity the gate. 
and the reasoning leading this equation indicate that orderly 
method determine experimentally the effect proximity the free surface 
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would investigate the flow phenomena with systematic variation 
and t/B. The writer has conducted series tests examine the mean flow 
pattern and the intensity turbulence (v’, the root-mean-square the velocity 
fluctuation) for the three conditions holding 2—namely, t/B 3.35, 
and The results are presented Figs. 32, 33, and The first value 
t/B was necessary for the formation free jump the vena contracta; the 
other values were chosen indicate the effects increased submergence. 


0.6 


Values 


Values 


enlightening consider the effect submergence from the standpoint 
the discharge characteristics the gate. Simultaneous consideration 
and the headwater depth, leads the expression: 


whenever appears any functional expression, Eq. 84, may re- 
placed Hence, the variation and h/B various degrees 
submergence should noted. dimensionless expression for the experi- 
mental determination this variation 


Experimental 
The discharge characteristics sharp-edged sluice gate, determined 
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experimentally the writer, are shown the solid curves Fig. 35. The 
center-line curves, lines constant were computed according Eq. 85. 
The three solid points represent values which the patterns mean flow 
and turbulence were investigated, and serve indicate clearly the change 
degree submergence. 

The experimental trend the curves, versus h/B, Fig. was verified 
approximate theoretical approach, utilizing the authors’ assumption 
hydrostatic pressure distribution throughout the flow plus the assumptions 
(see Fig. 31): energy loss between sections and uniform velocity dis- 
tribution sections and and, the vena contracta (section 2), zero 
velocity the roller and uniform velocity the jet. The energy relation 
between sections and and the momentum relation between sections and 
allow the following equations written: 


Values 


given the resulting forms Eq. for and for free discharge 
(assuming 0.6) are plotted Fig. the dotted curves. The dashed 
curve represents the results applying, for the common formula: 


which really the limiting form Eq. the depression the water surface 
decreases and approaches 

The writer’s data were obtained the laboratory the Iowa Institute 
Hydraulic Research (State University Iowa City), using sharp- 
edged sluice gate 1-ft flume which received water from constant-head 
tank. Velocities were determined pitot tubes the Prandtl and cylindrical 
type, the latter being used regions inclined flow. The fluctuations 
velocity due turbulence were determined turbulence (capable 
indicating only the logitudinal component), especially constructed the 
Iowa Institute for such measurements water. The results investigating 
the internal flow phenomena are shown Figs. 32, 33, and 34. The basic 
data these diagrams are the surface profiles, the velocity distributions, and 
Development Turbulence Pitot for Use Water,” Appel, thesis presented 


the State University Iowa, City, Iowa, June, 1949, partial fulfilment the requirements 
for the degree Master Science. 


Evaluation Several Experimental Methods Measuring Turbulence Flowing Water,” 
Hubbard and Norton, Report No. rept. prepared the Iowa Inst. Hydr. Research 
for the Office Naval Research under contract February, 1949. 


and 
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the distributions the intensity turbulence. All other data are computed 
from them. compared Figs. and 21, Figs. 32, 33, and indicate 
vividly the effect the proximity the free surface. 

direct result limiting the vertical extent the fluid, three basic 
changes occur the flow phenomena: 


(1) The discharge volume flux past successive sections perpendicular 
the jet, instead increasing continually according Eqs. and 44, must 
constant and equal the discharge the gate—that is, 


which the depth any section; 

(2) direct consequence item (1), fluid undergoing the induced 
motion because the accelerative action the jet must replenished 
return flow both and near the surface, forming what termed the 
and 

(3) The vertical expansion the jet does not continue indefinitely but 
limited the intersection the uppermost streamline the jet with the free 
surface. Inasmuch the free surface line constant pressure, this 
intersection must right angles, necessitating curvilinear pattern 
jet diffusion instead the linear pattern indicated Figs. and 21. The 
boundary the indicated the authors, corresponds 
closely the lines zero horizontal velocity Figs. 33, 34, and 35. 


The depression and unevenness the surface, due the gravitational 
effect, make invalid the reasoning leading they lead the 
constancy the sum pressure forces any vertical section and the flux 
the longitudinal component momentum the same section. Again, 
using the authors’ assumption hydrostatic pressure distribution (which 
the only one the three original assumptions that could reasonably apply 
flow conditions with finite tailwater depth) and letting represent the in- 
dicated sum, the following relation may written: 


which the specific weight water. The subscript Eq. applies 
uniform tailwater conditions. 

Actual values the ratios presented Eqs. and have been computed 
from the experimental data and appear above the mean flow patterns. The 
constant value unity indicated; the variations are doubt due part 
air entrainment and accelerative forces within the roller. 

With the complex flow pattern exhibited downstream from submerged 
sluice gate, one could not expect relations the nature Eqs. and 
represent energy flux any section. Instead, the writer has shown variation 


1r 
ti 
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total head, plotted above the mean flow patterns, computed from the 
experimental data follows: 


ce 


Here the authors’ assumption hydrostatic pressure distribution again 
utilized. 

any section the energy lost the mean flow, referred the energy 
the efflux section, the sum the kinetic energy turbulence and the 
energy expended heat over the preceding distance. This the relation 
indicated the authors (under the heading, the 


the kinetic energy the oncoming flow steadily [is] con- 
verted into kinetic energy turbulence, and the latter steadily 


This may expressed the equation: 


which the energy lost the mean flow; the kinetic energy 


turbulence; and the energy already expended heat. This relation 
represented schematically Fig. 31. 


The actual values were computed from the experimental data the 
following equation: 


69 


Eq. the assumption equal turbulence components the three coordi- 
nate directions has been made. The values computed were too small 
allow satisfactory representation indicated Fig. 31, but necessitated 
expanded scale plotted Figs. 32, 33, and discovery that the 
energy the turbulence, which involved great dissipation the mean 
flow energy, the order only great the total energy lost was one 
the most interesting results the writer’s experiments. 


allow rapid and vivid comparison, (the rate which energy 


lost the mean flow) was obtained graphically the slope the total head 
line and plotted the same set coordinates Because the smallness 
this rate represents essentially the rate that energy converted into 
heat. the scale the smallest turbulent eddies were constant throughout 


the flow, and should attain maximum values identical 


Since the writer had way determining the scale, such statement cannot 
made. The problem energy dissipation due turbulence will not 


completely solved until the discovery some method evaluating the size 
these smallest eddies. 
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1548 HENRY SUBMERGED JETS Discussions 


The writer has concluded that the method experimental attack 
applied the paper great value investigating any type submerged 
jet. only behooves the investigator recognize such factors the dissimi- 
larities resulting from boundary conditions, and modify the approach 
accordingly. 
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part the data for the mean flow patterns for 2.00, t/B and 
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COMING MEETINGS 


BOARD DIRECTION MEETING 


January 16-17, 1950. 
Quarterly Meeting will held New York, 
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HOTEL COMMODORE, NEW YORK, 
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The Reading Room the Society open from 9:00 5:00 Mondays 


through Fridays, inclusive. closed all day Saturdays, Sundays, and 
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Members, particularly those from out town, are cordially invited use this room 
their visits New York, have their mail addressed there, and utilize 
place for meeting others. There ample file current periodicals and 
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